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1.0  INTRODUCTION 

This  report  prepared  by  the  Link  Flight  simulation  Division  of  the  Singer 
Company  is  in  compliance  with  Contract  N61339-83-0068  granted  by  the  United 
States  Navy,  Naval  Training  Equipment  Center  (NTEC)  in  Orlando,  Florida,  to 
explore  the  implications  ot  Doppler  beam  sharpened  radar  and  synthetic 
aperture  radar  in  training  requirements. 

During  the  next  decade,  virtually  all  airborne  radars  will  have  Doppler  beam 
sharpening/synthetic  aperture  capabilities.  The  new  radars  ditter 
significantly  trom  present  real  beam  radar  and  perhaps  the  most  striking 
difference  is  in  the  respective  resolution  capabilities.  Rather  than  the  100 
to  300  toot  resolution  available  on  real  beam  radar,  synthetic  aperture  radar 
resolution  capabi lit ies  are  10  teet  or  less.  Within  the  next  tour  years,  it 
is  anticipated  that  three-toot  resolution  will  be  attainable. 

The  impact  that  accompanies  introduction  ot  the  new  types  ot  radar  is  not 
limited  to  resolution  factors  only,  dramatic  as  that  improvement  is.  In 
addition,  synthetic  aperture  radar  (SAR)  uses  principles  that  ditter  from 
current  real  beam  sets  and  consequently  the  radar  images  exhibit  different 
characteristics,  in  the  following  report,  these  differences  will  be  explored, 
anomalies  will  be  described,  and  the  impact  ot  the  differences  on  simulation 
will  be  discussed.  We  call  these  image  differences  anomalies  because  a  real 
beam  radar  s*"Ueni  and  a  synthetic  aperture  radar  system  would  generate  Images 
ot  dittering  appearance  even  it  both  systems  had  the  same  resolution. 

Apart  trom  the  resolution  improvement  and  the  anomalies,  the  new  radar  systems 
ditter  significantly  trom  present  ones  in  implementation  methodology.  While 


1.  There  is  no  universally  accepted  distinction  between  Doppler  beam  sharpened 
radar  and  synthetic  aperture  radar.  In  this  report  the  two  technologies  will 
be  used  interchangeably. 
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current  real  beam  radars  make  use  ot  some  digital  processing  techniques,  the 
new  generation  ot  radars  are  controlled  by  very  powertul  computers  with 
capabilities  tor  executing  several  million  instructions  per  second.  These 
advanced  computers  control  the  digital  signal  processing  devices  used  tor 
the  implementation  ot  synthetic  aperture  calculations  and  also  control  all 
modes  ot  the  radar.  The  resident  program  in  such  a  general  purpose  computer 
may  typically  consist  ot  several  hundred  thousand  instructions. 

Another  important  ditterence  between  real  beam  and  synthetic  aperture  radars 
are  the  failure  modes.  Aircrews  must  be  trained  to  recognize  these  tailure 
modes  and  make  appropriate  adjustments  to  restore  radar  system  tunctioning. 


The  introduction  ot  synthetic  aperture  radar  will  have  yet  another  major 
impact:  new  air-  to  ground  tactics  will  be  needed  because  the  resolution  ot 

SAR  images  in  the  vicinity  ot  the  t light  vector  degenerates  to  that  ot  real 
beam  radars.  Current  tactics  do  not  take  tull  advantage  ot  some  ot  the 
capabilities  ot  SAR  because  ot  this  "notch"  in  resolution.  There  is  no  doubt 
that  new  military  tactics  will  be  developed  which  use  this  high  resolution 
capability  to  increase  mission  effectiveness. 

In  presenting  the  current  SAR  Simulation  Study,  a  beginning  short  section  is 
devoted  to  current  radar  training  on  the  F-18  aircratt,  followed  by  a 
discussion  ot  SAR  theory,  radar  simulation,  and  relevant  discussions  ot 

anomalies  and  malt  unctions.  The  report  includes  photographs  ot  actual  SAR 

images  to  illustrate  various  ettects  and  characteristics.  These  images  are 
reproduced  by  courtesy  ot  Hughes  Aircratt  and  were  made  on  the  F-15  radar 

system  manufactured  by  Hughes  Aircratt.  The  report  also  contains  SAR  images 
made  by  the  Link  Emulation  Laboratory  from  an  enhanced  Detense  Mapping  Agency 
(DMA)  map  ot  the  China  Lake  area. 

As  will  be  explained  in  this  report,  radars  relying  on  the  Doppler  effect  have 
limitations  in  resolution  along,  and  in  the  vicinity  ot ,  the  flight  vector. 
This  area,  designated  as  the  "notch,"  is  typically  a  tew  degrees  in  width 
where  Doppler  Images  can  be  generated.  For  this  reason,  the  straight-in 

approach  toward  the  target  is  not  optimal  because  SAR  resolution  in  this  area 
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does  not  ditter  trom  that  ot  a  real  beam  radar.  Neither  does  the  straight- in 
approach  make  use  ot  one  ot  the  unique  teatures  ot  the  F/A- 18  radar;  that  is, 
rotating  the  display  in  such  a  manner  that  the  display  area  is  always 
correctly  oriented  with  respect  to  the  aircratt  heading.  Because  ot  these 
tacts,  it  is  expected  that  tactics  developed  tor  the  aircratt  will  be 
significantly  different  from  those  in  present  use. 

Because  synthetic  aperture  radar  simulators  are  likely  to  be  substantially 
ditterent  trom  present  ones.  Link  recommends  that  NTEC  procure  a  research 
synthetic  aperture  radar  simulator.  The  procurement  ot  such  a  simulator  would 
allow  the  U.S.  Government  to  accurately  determine  and  specify  requirements  ot 
future  simulators  and  also  assist  in  development  ot  new  tactics  tor  use  ot  the 
high  resolution  radar  systems.  The  research  simulator  should  be  highly 
programmable  to  allow  the  government  to  easily  perform  experiments  in  the 
selection  ot  special  effects ,  modes,  and  types  of  radar. 
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2.0  HIGH  RESOLUTION  RADAR  TRAINING 

At  present  the  F/A-18  is  the  only  operational  aircratt  with  high  resolution 
Doppler  beam  sharpening/synthetic  aperture  radar  capability.  To  determine  the 
status  ot  radar  training  on  this  aircratt.  Link  personnel  visited  both  LeMoore 
Naval  Air  Station  and  El  Toro  Marine  Corps  Air  Station  in  Calitornia.  LeMoore 
NAS  is  the  home  base  ot  an  F/A-18  Replacement  Air  Group  (RAG)  and  El  Toro  MCAS 
is  the  base  tor  an  F/A-18  Marine  Corps  squadron.  In  addition,  Link  has  met 
with  personnel  trom  Hughes  Aircratt  and  McDonnell  Douglas  Corporation  and 

discussed  the  technical  aspects  ot  radar  training  requirements  tor  this 
aircrat  t . 

At  present  the  F/A  18  squadrons  are  in  the  early  stages  ot  combat  readiness; 
their  principle  missions  are  detense  ot  the  tleet  rather  than  close  air 
support.  As  a  result,  training  is  primarily  concentrated  on  interceptor 
activities.  The  existing  F/A-18  simulator  has  radar  simulation  capabilities 
but  only  tor  air-to-air  engagements;  no  high  resolution  radar  landmass 

simulation  currently  exists.  To  train  pilots  in  the  use  ot  F/A-18  ground 
mapping  radar,  there  is  a  ground  school  program  tollowed  by  two  training 

tlights  to  teach  students  use  ot  the  radar  in  navigation  and  weapon  delivery. 

Considering  the  tact  that  most  students  have  had  no  previous  experience  with 

SAR  radar,  this  limited  amount  ot  training  is  insufficient  tor  pilots  to 

achieve  proticiency  in  the  numerous  ground  mapping  modes  and  capabilities  ot 
the  F/A  18  radar. 

It  appears  at  the  present  time  that  the  primary  training  requirement  tor 

pilots  is  tor  interception  modes,  with  a  secondary  goal  being  ground  support 
under  visual  conditions,  and  only  a  tertiary  requirement  tor  pilot  training  in 
close  air  support  under  limited  visibility  conditions.  Because  the  F/A-18 
aircratt  is  presently  only  a  single-seat  contiguration,  Navy  personnel  teel 
that  close  air  support  under  limited  visibility  conditions  is  teasible  only 
when  coordinates  oi  the  target  are  accurately  known.  Because  ot  the  task 
loading  on  the  pilot,  a  signit leant  amount  ot  radar  interpretation  during 

t light  is  not  teasible  and  hence  the  target  location  must  be  known  ahead  ot 


NAVTRAEQUIPCBNT  83-C-0068 


Because  oi  the  relatively  low  priority  ot  the  close  air  support  mission, 
tactics  tor  this  task  have  not  yet  been  tully  developed.  Squadron  VX-5  at 
China  Lake  is  now  tasked  to  develop  these  tactics.  It  is  clear,  however,  that 
these  tactics  must  ditter  trom  those  ot  current  aircralt  it  the  tull 
capabilities  ot  the  radar  are  to  be  used  to  support  the  attack  mission. 
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3.0  PRINCIPLES  OP  IK)PPI.RR  BEAM  SHARPENING  AND  SAR 

Numerous  articles  and  books  have  been  written  expressly  on  the  subject  ot 
synthetic  aperture  radar  (SAR)  and  collectively  they  cover  many  tacets  ot  this 
emerging  field  ot  radar  technology.  In  a  search  ot  this  literature,  however, 
no  single  article  would  lead  a  person  not  in  the  field  through  the  theoretical 
aspects  ot  SAR.  SAR  theory  is  founded  on  mathematical  concepts  and 
consequently  the  design  ot  a  SAR  system  requires  rigorous  analysis  and 
understanding  ot  fundamental  mathematics.  When  highly  advanced  mathematics 
are  required,  a  solution  is  available  from  accepted  and  traditional  sources. 
For  simulation  purposes,  comprehension  of  SAR  technology  is  essential  tor 
discriminating  selection  ot  important  factors  and  characteristics,  and  also 
tor  recognizing  traits  that  are  irrelevant  in  radar  simulation. 

In  this  SAR  study,  theory  development  progresses  by  using  basic  mathematics 
and  physics  laws  and  supplemental  sources,  as  required.  Reference  sources  are 
listed  and  identified  in  Appendix  A.  Subjects  in  the  study  include  Doppler 
theory,  synthetic  aperture  radar  theory,  equations,  simulation  algorithms,  and 
radar  simulation  algorithms.  A  list  of  mathematical  signs,  definitions,  and 
symbols  is  included  in  Appendix  B. 

3 . 1  DOPPLER  EFFECT 

Doppler  shitt,  the  change  in  frequency  of  a  wave  reaching  an  observer  or 
system,  is  the  basic  concept  tor  synthetic  aperture  radar  and  a  logical 
starting  point  tor  the  SAR  study.  Doppler  equations  are  derived  trom  two 
postulations.  The  first  approach  is  classically  derived  and  the  second  is 
distinguished  by  use  ot  the  Doppler  shitt  in  SAR  applications.  In  both 
derivations,  an  aircraft  will  be  assumed  to  be  heading  directly  toward  a 
reflecting  surface. 

In  the  classical  approach,  a  vehicle  transmitting  at  a  frequency  F  and 

T 

moving  in  the  direction  ot  transmission  will  produce  F  cycles  in  a  length 

IC-V  J  where  C  is  the  speed  ot  light  and  V  the  velocity  ot  the 
AC  AC 

aircraft.  The  wavelength  ot  this  signal  will  be  l(C-V  )/F  J.  To  an 

AC  T 

observer  on  the  ground,  the  received  frequency  will  be  IC/(C-V  )/F  J  or 

AC  T 
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IT"  0T.n'  - 
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(OF  /(c-V  )J.  This  trequency  is  that  ot  a  retlected  signal  trom  a 

T  AC 

vehicle  approaching  a  retlecting  surtace.  Assume  that  the  transmitted  signal 

is  ahead  ot  the  vehicle.  The  velocity  ot  the  passing  signal  is  C+VAC*  The 

number  ot  wavelengths  passing  per  second  is  l(C+V  )/CJ,  hence  the  received 

AC 

trequency  will  be  IF^IC+V^J/CJ .  It  we  summarize  the  above  equations  and 
combine  them,  the  equation  tor  the  transmitter  moving  toward  a  stationary 
observer  is: 


1  **•' 
’">"!*•***• 

■/V.v'v 


F»  =• 


(3.1.1) 


Where 


F  =  transmitted  trequency 

Fr  =  received  trequency 

C  =  velocity  ot  light 

V  =  aircratt  velocity 

AC 

It  the  transmitter  is  stationary  and  the  observer  is  moving  toward  it, 
the  equation  is: 


Ft\C+Vac\ 


is* 


*  >  .*•  ‘ 


(3.1.2) 


The  equation  tor  a  Doppler  trequency  it  the  transmitter  is  moving  and  the 
signal  is  retlected  back  to  the  transmitter  is  obtained  by  substituting  the 
received  trequency  ot  Equation  3.1.1  tor  the  transmitted  trequency  in  Equation 
3.1.2. 


*  s  *  *■  *  .  *, 


(C+V^) 

(C-  VAC) 


(3.1.3) 


The  Doppler  shilt  is  obtained  by  subtracting  the  transmitted  trequency  trom 
Equation  3.1.3: 


Fps  -2  Vac  Ft!(C-  V„c) 


(3.1.4) 


/.*\V 
O  S’  */.« 
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Where 

F  =  Doppler  trequency 

Note  Equations  3.1.1  and  3.1.2  give  a  slightly  ditterent  trequency  shitt. 

The  SAR  approach  does  not  ditterentiate  between  the  outward  path  and  the 
return  path  ot  the  signal.  Consider  the  signal  ot  the  transmitter  to  be 
represented  by  the  signal: 

Sr  »  />  /]  (3.1.5) 

Where 

exp  =  exponential  to  the  base  e 

j  =  square  root  ot  -1 

«  =  3.14159 

t  =  time 

S  =  transmitted  signal 

T 

The  range  trom  the  transmitter  to  the  ret  lector  can  be  given  by: 

R  —A, |  +  Rr  t  (3.1.6) 


Where 


R  =  range 

R  .  =  range  at  time  tl 
tl 

R  =  range  rate 
R 


The  elapsed  time  tor  the  signal  to  reach  the  retlector  and  return  is: 
(/?,.  +  Rpt) 


fa 


=  2  R/C 


r. 


(3.1.7) 
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Where 


t  =  delta  time 
D 


In  Equation  3.1.5  substitute  t+tQ  tor  t  to  give: 


S„  *=  tip 


.  „  „  4  .  }  *  *  FT  RR  t  jA-n-FrRn 

>**-Fr*  + - c - + - c - 


(3.1.8) 


Where 


S  =  received  signal 
R 


The  term  j'4'vF^R^/C  is  a  constant  phase  shltt  and  can  be  dropped 
out  to  give: 


4*jr  ‘F't'Rr) 

ST=txp  jt(2x  FT  + - - - 


(3.1.9) 


Equation  3.1.9  can  be  divided  into  two  parts  by  using  Equation  3.1.10: 


txp{jz)  =c0s(2)  +  j-sin(jr) 


(3.1.10) 


Where 


cos  =  cosine 
sin  =  sine 
X  =  variable 
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AifFj-Rfft  .  .  4 

SH  =coi  2  k  Fj  t  + - - -  +  ;sin  2  k  Ff't  + - - - 

c  c 


(3.1.11) 


Subtracting  the  transmitted  signal  Irom  Equation  3.1.11  and  substituting  the 
velocity  ot  light  divided  by  wavelength  tor  transmitted  trequency  gives: 


4  K  -Rp  t  i-K  -Rfi  t 

Sfi  —  cos - - -  +  j  sin  - - - 

A  A 


(3.1.12) 


Where 


X  =  wavelength  ot  the  signal 


The  received  signal  is  in  units  ot  radius  per  second  and  range  rate  (R  )  is 

R 

the  aircratt  velocity.  The  signal  can  be  expressed  in  cycles  per  second  by 
removing  2*ir.  This  simplit ication  is  given: 


c  2  VAC  t  _  .  .  2  VAC  t 

Skf=C03 — - —  +jsinl — - — 


(3.1.13) 


Where 


S  =  received  signal  in  trequency 
Kr 


When  the  in  phase  and  quadrature  phase  are  sampled  at  twice  the  trequency  ot 
the  received  signal  trequency  or  higher,  the  complete  description  ot  the 
retlected  signal  is  known.  The  amplitude  ot  the  retlected  signal  is  given: 


—  V  Sip  +  Sn 


(3.1.14) 


Where 


S  =  received  signal  amplitude 
AMT 

Sip  =  in  phase  component  ot  the  received  signal 
3^  =  quadrature  component  ot  the  received  signal 
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The  phase  shitt  ot  the  reflected  signal  is  given: 


P 


(3.1.15) 


Where 

tan  =  tangent 

In  operational  SAR*  the  retlected  signal  is  converted  to  digital  signals  by 
sampling  the  retlected  signal.  This  sampling  is  normally  perlormed  at  a  rate 
that  is  higher  than  the  Nyquist  sampling  rate.  The  number  ot  levels  in  the 
analog-to-digital  converter  determines  the  dynamic  range  ot  the  received 
signal.  When  the  digital-to-analog  converter  saturates,  intormation  is  lost 
and  ghost  signals  can  be  generated.  Dynamic  range  and  saturation  tactors  are 
discussed  in  other  sections  ot  this  study. 

In  our  reviews,  it  has  been  shown  that  by  sanjpling  the  in  phase  and  quadrature 
phase  ot  a  retlected  signal  at  a  sampling  rate  ot  at  least  twice  the  retlected 
signal  frequency,  all  characteristics  ot  an  image  point  can  be  recorded. 
Recording  of  the  reflected  signal  to  allow  signal  processing  to  be  performed 
in  non-real  time  is  the  basis  of  modern  synthetic  aperture  radar. 

3.2  SAR  THEORY 

The  theory  of  Doppler  shitt  has  been  briefly  introduced.  The  theory  ot 
synthetic  aperture  radar  is  developed  using  the  Doppler  equation.  Rather  than 
straight  ahead,  the  reflector  will  be  90  degrees  from  the  velocity  vector. 
Results  when  the  target  is  not  90  degrees  from  the  velocity  vector  will  be 
developed  in  later  text.  Figure  3.2-1  is  a  diagram  used  to  develop  SAR 
formulations. 

The  line  D  (distance  from  point  PT)  is  perpendicular  to  the  velocity 
rPT 

vector  and  connects  the  velocity  vector  to  point  P  (X  ,Y,  This 

T  (n)  (m) 

line  is  the  reference  line  of  the  SAR  equation.  The  signal  that  is  received 
after  being  reflected  from  the  ground  is  given  by: 
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S* 


^Sj-Tusr  [Z(n)>RsHm))'e*P 


(3.2.1) 


Where 


T  =  target  retlectance  coetticient 

REF 

X  =  azimuth  coordinate 

(n) 

R  =  range  coordinate 

SUm) 

F  =  frequency 

R  -  slant  range 

SL 


The  slant  range  (R  )  can  be  derived  from  the  distance  (D„  )  and  the 

SL  FPT 

distance  (Y-Y  ): 

(m) 


Rsi  —  \f  [Dfpt]  +  (*"  -  ^(*))2  (3.2.2) 


Where 


D  =  distance  trom  point  PT 
FPT 

Y  =  reterence  along  the  flight  path 

Y  *  point  along  the  flight  path 
(m) 


Since  the  term  (Y-Y  )  is  much  less  than  D  ,  the  binomial  expansion  can 

(m)  FPT 

be  used  to  get  a  simpler  expression: 


Rsl  —Dfpt  +  - — -  -  (3.2.3) 

*  'VrPT 


Equation  3.2.4  is  obtained  by  substituting  Equation  3.2.3  into  Equation  3.2.1: 

2  Dr„/C  -  |K-  Y{m)]2JC  DrPT) J  (3.2.4) 


3-8 
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&  jr^ 

*•  -  wzl 


Collecting  the  terms  in  Equation  3.2.4  gives: 

£jt(a.a)  ^(*(,),/?5i( *)) j-  ezp^j-2*  -F-(  -  4  X  Dfrr  F/C  - 

2nF(Y-  K(b))3  !C-D„t] 


(3.2.5) 


Substituting  the  velocity  ot  light  divided  by  wavelength  tor  trequency  and 
separating  the  exponent  gives: 


eip^-  j ■(■»  D^.r/xj  erp[-  j  (2  x  •(  Y-  r(B))  A  Dfpr\ 


(3.2.6) 


The  tirst  exponent  is  the  transmitted  signal  and  contains  no  intormation.  The 
second  exponent  is  a  constant  phase  shitt  and  contains  no  intormation.  By 
demodulating  the  signal,  the  tirst  and  second  exponent  can  be  removed: 


—  st  tret  j-2  n  -[Y  -  l<m))  A'^V/rj 


(3.2.7) 


£ 


AU 


It  the  tactor  2  is  removed.  Equation  3.2.7  would  be  the  equation  ot  a  positive 

cylindrical  lens  centered  at  y,  .  Because  the  SAR  antenna  transmits  and 

(m) 

receives  the  signal,  there  is  a  two-way  path  tor  phase  shitt  occurrence; 
whereas,  a  lens  has  an  external  source  ot  illumination  and  there  is  only  a 
one-way  path  tor  phase  shitt  occurrence.  Equation  3.2.7  can  be  simplitied  to 
in  phase  and  quadrature  components. 

This  simplit icat ion  would  be  a  SAR  tunction.  For  this  treatise,  the  equation 
as  given  is  more  meaningtul. 

2-9 
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In  Equation  3.2.7  each  signal  (S  )  is  composed  ot  the  sum  ot  all  the 

R(n) 

reflected  energy  trom  the  points  on  the  arc.  When  the  sample  ot  the  arc  trom 
points  U1  to  U2  is  made,  the  retlectance  ot  the  arc  becomes  known. 
Intormation  is  in  a  coded  format,  and  to  be  usetul  it  must  be  decoded.  In 
optics,  this  coded  intormation  is  reterred  to  as  a  holograph.  The  tirst  SAR 
processors  made  use  ot  holographic  principles  to  decode  SAR  images.  To  obtain 
SAR  images,  azimuth  and  slant  range  must  be  scanned. 

To  store  SAR  image  Intormation,  two  memories  are  required.  one  memory  is 
used  to  record  the  in-phase  (cos)  signals  and  the  other  memory  is  used  tor 
recording  the  quadrature  signal  (j*sin).  The  memories  can  be  considered  to 
be  two-dimensional.  Loading  takes  place  along  the  slant  range,  a  single 
scanline  at  a  time.  When  data  is  decoded,  access  is  at  a  constant  slant 
range.  This  type  ot  memory  is  otten  called  a  corner  memory.  This  theoretical 
account  has  identified  none  ot  the  second  order  effects  which  later  will  be 
defined  individually.  The  report  now  equates  synthetic  aperture  radar  wi'h 
the  nature  and  properties  ot  optics. 

A  valid  reason  tor  analyzing  SAR  and  its  relationship  to  optics  is  that  a  lens 
is  a  commonly  known  device  familiar  to  most  everyone.  Figure  3.2-2  is  a 
diagram  ot  a  1:1  cylindrical  lens.  The  cylindrical  lens  collects  rays  trom 
point  0  trom  the  object  and  redirects  them  to  form  point  P  at  the  image 
plane.  In  the  center  ot  the  lens,  rays  trom  a  number  ot  points  in  the  object 
plane  meet  at  a  single  point  in  plane  L  of  the  lens.  This  happens  along  all 
points  in  plane  L.  A  SAR  system  duplicates  the  retraction  of  the  cylindrical 
lens  trom  the  object  to  plane  L.  SAR  stores  intormation  as  it  appears  at 
plane  L  and  processes  it  later;  a  lens  system  performs  real-time  processing. 

Some  ot  the  tirst  SAR  processors  used  an  optical  system  tor  processing  SAR 
images;  this  type  ot  processing  is  still  continued.  A  good  discussion  ot 
processing  SAR  data  by  optical  methods  can  be  found  in  "Introduction  to 
Fourier  Optics"  by  j.w.  Goodman  (page  184).  Modern  SAR  systems  in  aircrait 
process  by  high  speed  signal  processing.  Some  ot  the  problems  evident  in  SAR 
can  be  observed  by  using  a  lens.  By  blocking  a  part  ot  the  lens,  the  image 
dims  and  is  less  sharp,  causing  lens  aberrations.  In  SAR,  acceleration  errors 
do  the  same  thing  to  a  SAR  image.  When  a  prism  is  put  between  the  object  and 
the  lens,  the  image  is  shitted.  In  SAR  a  velocity  error  will  cause  the  SAR 


LENS 


DIAGRAM  (CYLINDRICAL  LENS) 
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image  to  shitt  in  the  same  way  as  the  prism  in  the  optical  system.  From  this 
explanation,  it  can  be  seen  that  there  is  nothing  mysterious  about  SAR. 


The  next  step  in  understanding  synthetic  aperture  radar  is  more  complex.  The 
book  "Introduction  to  Statistical  Optics,"  by  Edward  L.  O’Neill  has  a  chapter 
on  Dittraction  Theory  ot  Image  Formation  (page  70).  In  this  chapter,  tour 
equations  are  developed  tor  a  cylindrical  lens  and  the  example  ot  a 
cylindrical  lens  can  be  related  to  the  equations  ot  SAR  tor  better 
understanding.  The  tirst  equation  is  the  complex  amplitude  distribution  in 
the  exit  pupil.  The  intormation  trom  this  equation  reflects  what  is  stored  in 
SAR  memory.  This  equation  is: 

F(P)  —  0  c  £  (inride  the  aperoture}  ^  2  8) 

=  0,  p  c}.  £  (outride  the  aperoture) 


P 


k 


2n 

X 


Where 

k 

A 

F(B) 

i 

K 


=  wave  number 
=  lens  aberration 
=  aperture  function 
=  square  root  ot  -1 
=  selectable  constant 


(physics  usage) 


The  second  equation  is  the  complex  amplitude  distribution  in  the  point  (line) 
image : 


“(O 


00 

=  AJ  F(0)e-"'dp 

-  OO 


(3.2.9) 


Where 

u 

A 


u 


(X) 


=  dimension  in  the  aperture 
=  amplitude 

=  complex  amplitude  distribution  of  the  point  image 
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Equation  3.2.9  describes  the  distribution  ot  light  trom  a  point  source  at  ter 
going  through  a  lens.  It  the  lens  were  pertect,  a  point  would  remain  a 
point.  This  equation  is  very  important  and  it  will  be  used  to  derive  the  SAR 
resolution  equation. 


The  third  equation  is  the  intensity  spread  tunction: 


(3.2.10) 


Where 

S(x)  =  intensity  spread  tunction 

Equation  3.2.10  is  a  description  ot  the  energy  distribution  ot  a  point  source 
alter  going  through  a  lens. 

The  tourth  equation  is  the  normalized  transter  tunction. 

3.2.11  and  3.2.12: 

/»(*)*-'  dx 

/  V  “  00 

r(Wj)  - - - - 

/  »(x)rf* 

-  oo 

Where 

t  =  normalized  transter  tunction 
(ux) 


These  equations  are 


(3.2.11) 


OO 

/  F(0)F’i0  -  u,)df) 

-  00 

00 

/  \F(P)\2d0 

-  00 


(3.2.12) 


Where 

6  =  K*»i/R 

w  =2**  times  trequency 
R  =  distance  trom  aperture  to  image  plane 
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The  derivation  ot  the  character istics  ot  a  cylindrical  lens  can  be  tound  on 
page  79  ot  "Introduction  to  Statistical  Optics."  The  response  distribution  ot 
a  point  Is  given: 


u(x)  —2A 


6in(0o*) 


(3.2.13) 


Where 


u  =  response  distribution  ot  a  point 
(X) 

a  =  one-halt  ot  the  aperture  width 

A  *  amplitude 


The  normalized  spread  t unction  is: 


*(j)  = 


6in(0ox 


(3.2.14) 


The  transter  tunction  is  given: 


'<"■>=  '-5K- 


I",  I  <  w. 


(3.2.15) 


'v  • 


=  0 


|w,  |  >  20o 


Figure  3.2-3  is  a  plot  ot  these  tunctions.  The  equations  ot  the  lens  can  now 
be  used  to  derive  the  equations  ot  SAR. 


The  resolution  will  be  derived  trom  Equations  3.2.13  and  3.2.15.  By 
inspection  it  can  be  seen  that  Equation  3.2.13  is  zero  at  «,  2%  ...n«  . 


«  _  71 

0o  2  *-o/R-\ 


(3.2.16) 


.-V 

l  -  * 


R\  _  R  \ 

2fl  Ai 


(3.2.17) 
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Where 

A^  =  aperture 


An  interpretation  ot  Equation  3.2.17  shows  that  a  point's  main  intensity  area 

is  spread  ovei  a  distance  ot  2RX/A.  In  Equation  3.2.15  let  r(w  )  be 

x 

zero  and  solve  tor  u  . 

x 


2-ft  Al 

fix 


(3.2.18) 


The  wavelength  in  the  transfer  function  that  has  zero  amplitude  in  Equation 
3.2.18  is  indicated  in  Equation  3.2.19.  the  resolution  tor  a  lens.  It  has 
already  been  mentioned  that  range  would  be  divided  by  two  (the  two-way  path). 
Hence  the  resolution  tor  SAR  can  be  written  as: 


RES  =  — — 


(tor  a  lens)  (3.2.19) 


RES  — 


fiX 
2  Al 


(tor  SAR)  (3.2.20) 


Where 

RES  =  resolution  in  length 

This  resolution  occurs  when  the  modulation  transfer  function  has  zero 
modulation.  Two  pixels  are  required  tor  each  cycle  by  the  Nyquist  criteria. 

Equation  3.2.20  is  the  resolution  equation  ot  a  SAR  that  is  looking  at  the 
ground  90  degrees  from  the  velocity  vector  and  "A"  is  the  data  collection 
distance  along  the  flight  path.  Figure  3.2-4  is  a  diagram  indicating  what 
happens  when  data  tor  SAR  processing  is  not  collected  at  90  degrees  trom  the 
velocity  vector.  "A"  is  the  data  collection  distance  but  the  effective  data 
collection  distance  along  the  velocity  vector  is  A*sin($).  SAR  processing 
is  most  etticient  at  90  degrees  and  is  zero  straight  ahead. 
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Equation  3.2.20  can  now  be  written  as: 


EES 


R  SL 

2£>csin(  <p) 


(3.2.21) 


Where 

=  distance  along  the  (light  path  where  data  is  collected 
$  *  squint  angle 

Equation  3.2.20  can  be  derived  in  other  ways.  Another  way  is  to  differentiate 
the  Doppler  frequency  with  respect  to  the  squint  angle: 

(3.2.22) 

(3.2.23) 

Equation  3.2.23  shows  that  SAR  processing  Is  most  efficient  at  90  degrees  from 
the  velocity  vector  and  the  resolution  Is  a  function  of  the  sine  of  the  squint 
angle. 

Information  and  equations  tor  SAR  simulation  requirements  are  now  stated. 
Equation  3.2.9  shows  that  the  designer  of  a  SAR  system  has  a  powerful  tool 
with  which  to  control  the  side  lobes  of  the  radar.  Note  that  response  of  a 
point  is  a  function  of  the  aperture  (F(0)).  When  the  samples  are  multiplied 
by  a  weighting  function,  the  response  of  a  point  can  be  modified.  Figure 
3.2-5  is  a  plot  ot  uniform  aperture  and  a  weighted  aperture.  By  accepting  a 
wider  main  lobe,  the  side  lobes  can  be  decreased. 


Fp  =Ft  cos(<s) 
dFR 

— —  =Ft  sid(<S) 


There  are  many  functions  that  can  be  used  to  modify  the  aperture.  One  such 
function  is  the  cosine  raised  to  the  k  power.  Figure  3.2-6  is  a  plot  of  the 
point  response  as  a  function  ot  k  in  which  the  function  is  cos* 
Modifying  the  aperture  with  the  function  cos  attenuates  the  higher 
frequencies  ot  the  normalized  transfer  function. 
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Since  the  aperture  function  is  being  modified  to  control  the  side  lobes,  there 
should  be  a  better  method  of  specifying  resolution.  One  way  to  do  this  is  to 
take  into  account  the  3  dB  bandwidth.  Equation  3.2.21  can  be  rewritten  to 
accomplish  this,  as  indicated  in: 


RES 


Rsl  ^  hr 
2  Dc  nia(?>) 


(3.2.24) 


Where 

K  =  aperture  modification  factor 

F 

4 

When  the  aperture  function  is  cos  ,  then  k  would  be  1.94.  This  value  is 

r 

twice  as  much  as  the  limit  of  the  normalized  transfer  function.  Table  3.2-1 
is  a  comparison  of  real  aperture  radar,  SAR,  and  a  lens  function. 

3.3  ISODOPPLER  MAPPING 

Theoretically,  SAR  was  developed  by  considering  only  one  arc.  This  is  correct 
for  SAR  processing  tor  each  arc  contains  all  data  necessary  to  recover  the 
azimuth  information.  Unfortunately,  arcs  are  not  linear  with  range.  Figure 
3.3-1  is  an  exaggerated  example  of  isodoppler  nonlinear  mapping.  This  figure 
shows  a  20  degree  sector  scan  that  is  centered  at  an  azimuth  of  30  degrees. 
At  mid-range  the  isodoppler  filters  were  set  to  match  the  sector  scan.  The 
shaded  area  is  mapped  with  constant  isodoppler  tilters.  This  phenomenon 
distinguishes  synthetic  aperture  radar  from  real  aperture  radar. 

To  understand  isodoppler  mapping,  an  error  equation  is  derived  and  examples 

are  given.  Let  K  -  cos(0  )»cos(¥  )  where  6  and  T  are 

R  R  R  R  R 

reference  angles  and  the  isodoppler  tilters  match  the  width  of  the  sector 

scan.  The  equation  tor  azimuth  error  mapping  is: 


©r>r  *=©*  -  cos"  *  f—  ,R  r  1 


(3.3.1) 


NAVTKAEQU1PCEK  83- C- 0066 


Z  *<  L/O  I/O  d 

*—  ^  'v*  LU  u_ 

la  cr  «—  z  o  x 

I/O  '-V.  >•  LU  — 1  UJ 

uj  ^  r  2  _ >  oo  — 

x  d  c-^  — «  x  z 

c~c?  —  cc  u  -  r  'J 

cr  lo  ui  c  -  3  z 

a.*-'  c  _j  o  -i  u  uj  z  « i 

Z  C"u  d  d  CL-  <__>  d  LA 

_ I  •— «  *  I  UJ  CD  a  Lu  lO 

cr  oo  i  cr  c  cd  uj  u.-  -» 

Z  *  —  >-  O  Cl  X  O  O 

Ucra-zcroozzoro 
— •  *  xOLuOoOdddcr 

i/,  c\ujZ>o*-ucr-(i. 


o  z 

— «  <r 

z  *—  u-  x 
C  «s  tr 

—  x  cc 
la:  a-  d  >- 

UJ  U  L  CD 


d  X>  Q_  Z  X 

Z  d  Q.  ■—  H- 

Z  X  d  _J  2D 

u.  lu  uj  Z  Q:  2 
x  —  *-dOO  — 

lu  z  o  o  x*  cr  o  ^ 

lA  d  Z  Z  LA  CD  Q_  d 


^  —  —  cr 

^  O  2T  *—  u_ 

Q_  D  u-  O 

>■  x  Z 

TD  *—  _j  —  21 

—  Z  d  <-A>  •— «  X  * — 

_  lu  o  o  uj  a. 

—  c  z  o  ^  o  a;  j  u.  u 

d  CNj  ►— *  O  LU  LU  d  O  O 

z  — >►-  u  u  a  ^  z  2 

cr  i  lu  uj  z  cr  x 

Lu  LU  *<  • — ■  — I  UJ  >•  n  t-  >-1  UJ  _ X 

►—zcrxdzcroz  o 

XO'^XLUOlUCTUJU.X  —  _J 
uz^ua:z>Q.uOuiij 


I/O 

x> 

1  < 

*— 

LO 

cr 

. _ / 

LU 

a 

* — 

cr 

LU 

X> 

>- 

L— 

X 

d 

» — 

z 

*— H 

cr 

►— * 

a 

z 

o 

♦— 

d 

- 1 

LU 

X 

(_■ 

L— 

X 

d 

O 

z 

Lu 

X 

; _ ' 

X 

d 

»— 

Lu 

o 

z 

c? 

►— 

cl 

cr 

*— > 

X 

—* 

O 

z 

T" 

c: 

wo 

lO 

d 

LU 

X 

d 

— 

lO 

o 

cr 

z 

> 

_J 

LU 

Lu1 

_J 

X 

I— « 

tA 

C- 

■X 

cr 

LA 

cr 

»— 

c. 

*~~l 

d 

Lu 

\j 

d 

o 

LU 

c — 1 

Cu 

X 

L  O 

— I 

UJ 

~7* 

cr 

>• 

•— « 

d 

c 

— 

■ — 

<- Ai 

Cr 

r_ 

Cl. 

X 

Cu 

X 

u_ 

NAVTRAEQU1PCENT  83-C0068 


Where 

6  =  azimuth  reterence  angle 

R 


k  =  cos  (6  )*cos(¥  ) 
R  R  R 


♦  =  depression  angle  reterence 


¥  =  depression  angle 


0  =  azimuth  Doppler  error 

DE 


Dittexeutiat ing  Equation  3.3.1  with  respect  to  depression  angle 
3.3.2: 


d&DE  kif  sin(V’) 


dtp  cos(V’) 


Vcos2(V>)  -  t/l 


Equation  3.3.2  can  then  be  simplltied: 

dOpe  _  tanfV’/?) 
dtp  tan(0fl) 


Figure  3.3-1  contirms  what  Equation  3.3.3  implies.  General 
occur  at  large  depression  angles  and  near  the  velocity  vector 
given: 

Example  One:  Isodoppler  Mapping 
Center  ot  scan  is  30  degrees 
Depression  angle  reterence  is  3  degrees 
Depression  angle  ottset  is  1  degree  and  3  degrees 


gives  equation 

(3.3.2) 


(3.3.3) 


mapping  errors 
Examples  are 
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1  degree  ottset  causes  -0.120688  degree  error  in  azimuth 
3  degree  ottset  causes  +0.243451  degree  error  in  azimuth 

Example  Two:  Isodoppler  Mapping 

Center  ot  scan  is  30  degrees 

Depression  angle  reterence  is  10  degrees 

Depression  angle  ottset  is  5  degrees  and  15  degrees 

5  degree  ottset  causes  -1.115607  degrees  error  in  azimuth 
15  degree  ottset  causes  +2.000800  degrees  error  in  azimuth 

In  the  design  ot  a  SAR  simulator,  it  is  necessary  to  make  a  determination  ot 
isodoppler  frequency  in  the  aircratt  SAR  system.  In  radar  simulators,  tor 
either  synthetic  or  real  aperture  systems,  aspect,  shadow,  and  slant  range 
processing  must  be  performed  along  scanlines.  This  is  an  important 
consideration  when  designing  tor  isodoppler  mapping.  The  tact  that  isodoppler 
mapping  is  nonlinear  with  range  requires  interrogation  ot  more  scanlines  in 
synthetic  aperture  radar  than  in  real  aperture  radar  because  some  ot  the  data 
generated  in  SAR  will  not  be  used.  Figure  3.3-2  is  an  example  ot  SAR 
simulation  overscan.  The  shaded  area  represents  data  that  will  not  be  used  in 
the  display. 

The  method  ot  data  display  is  determined  by  the  aircratt  SAR  system.  When  the 

aircratt  correctly  maps  data,  simulator  data  is  already  correctly  mapped. 

Thorn  is  a  problem  in  correctly  determining  what  data  is  in  the  shaded  area 

and  assuring  that  this  data  is  not  displayed.  It  the  aircratt  does  not 

correct  tor  isodoppler  mapping,  the  simulated  isodoppler  mapping  becomes  more 
1 

complex  .  The  nonlinear  isodoppler  scan  must  now  be  mapped  into  the  linear 
sector  scan.  Figure  3.3-3  is  a  diagram  ot  this  mapping. 


1.  F/A-18  radar  isodoppler  mapping  errors  are  corrected  when  data  is  written 

on  the  display.  This  means  that  the  display  is  a  distorted  sector  scan. 


Figure  3.3-2  ISODOPPLER  SCAN  RELATIVE  TO  SECTOR 
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From  a  theoretical  point  of  view,  there  is  no  reason  which  prevents  generation 
ot  a  frequency  (a  function  ot  range)  to  correct  tor  isodoppler  frequencies. 
This  generated  frequency  could  be  used  as  the  local  oscillator  frequency  to 
bring  Doppler  frequencies  to  SAK  processing  frequency.  It  the  aircraft  SAR 
system  performed  in  this  way,  simulator  isodoppler  mapping  would  no  longer  be 
a  design  problem.  At  the  current  time,  limited  data  on  isodoppler  mapping 
prevents  a  detailed  analysis  and  satisfactory  solution  to  design  ot  simulator 
isodoppler  mapping.  Items  of  consideration  in  designing  a  simulator  SAR 
isodoppler  module  include: 

1)  Performance  characteristics  of  the  aircraft  SAR  system 

2)  Accuracy  versus  training  value 

3 )  Cost 

3.4  SQUINT  ANGLE  DEPENDENCE 

Squint  angle,  the  angle  between  the  physical  axis  ot  the  antenna  and  the  axis 

ot  the  radiated  beam,  enters  into  many  SAR  calculations  and  is  a  good 

beginning  for  deriving  SAR  equations  for  simulation  purposes.  Figure  3.4-1  is 

a  diagram  ol  the  squint  angle  ($),  the  azimuth  angle  (8^),  the 

depression  angle  (T),  the  altitude  (A),  the  ground  range  (R  ),  and  the 

'  GR 

slant  range  (R  ).  Appendix  B  contains  a  definition  ot  the  above  mentioned 
SL 

terms.  Using  trigonometric  functions,  the  following  relationships  can  be 
derived  relative  to  Figure  3.4-1: 

X  =/?c-8in(ew)  (3.4.1) 


(3.4.2) 


«sl 


Rgh 

cor(  V’) 


(3.4.3) 


COS  (6yv) 


(3.4.4) 


**-  -  1 
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cos(V’)  — 


Vx^TWTa^ 


cos(e) 


Y 

y/X2  +  Y1  +  A2 


cos(0)  =  cos(9yV)cos(4' ) 


(3.4.5) 


(3.4.6) 


(3.4.7) 


Where 

0  =  azimuth  angle 

vv 

A  =  altitude 

R^  -  ground  range 

Equation  3.4.7  Is  derived  by  multiplying  Equation  3.4.5  and  Equation  3.4.6. 


Figure  3.4-2  is  derived  by  passing  a  plane  through  the  velocity  vector  (V  ) 

v 

and  the  slant  range  (R  ).  Velocity  vector  dimensions  are  length/time,  and 

SL. 

slant  range  dimension  Is  length.  Only  directions  ot  the  vectors  were  used  to 
derive  the  plane.  The  closing  velocity  Is  defined  as: 

V'c  *  Vvcos(<t>)  (3.4.8) 


Where 

V^  =  closing  velocity 

V  =  velocity  vector 

v 

in  Figure  3.4-2  v^  Is  the  closing  velocity  ot  the  point  (P)  and  the  alrcratt 
(AO.  When  P  is  straight  ahead  ot  the  flight  path,  the  closing  velocity 
(V^)  is  equal  to  the  velocity  vector;  when  the  point  is  abreast,  the  closing 
velocity  is  zero.  The  Doppler  frequency  and  the  closing  velocity  are  directly 
related  and  thus  the  solution  tor  the  closing  velocity  is  also  the  solution 
tor  the  Doppler  frequency. 
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In  synthetic  aperture  radar/  a  point  on  the  ground  Doppler  trequency  changes 
as  the  aircratt  tlies  along  the  data  collecting  path,  in  analyzing  SAR  tor 
simularion  purposes,  the  Doppler  trequency  ot  a  point  haltway  along  the  data 
gathering  path  is  called  the  Doppler  trequency  signature.  It  we  return  to 
Figure  3.4-1  and  the  associated  equations  and  combine  Equations  3.4.6  and 
3.4.8,  the  resulting  equation  is: 


Vc  = 


VyY 

v/A^  +  K8  +  A 2 


(3.4.9) 


Equation  3.4.9  can  be  normalized  by  setting  and  A  to  1. 


(3.4.10) 


Equation  3.4.10  can  be  written  in  the  torm: 


-A*  + 


i  r2  =i 


(3.4.11) 


Equation  3.4.11  is  the  equation  ot  a  hyperbola  symmetric  to  the  Y  axis.  By 

setting  V  to  dil Inrent  values,  a  tamily  ot  hyperbolas  can  be  generated, 
c 

When  Y  is  positive,  Doppler  frequencies  are  tor  a  closing  object  and  when  Y  is 
negative,  Doppler  trequencies  are  tor  an  opening  object.  Figure  3.4-3  is  a 
diagiar.)  ot  isodoppler  and  isorange  lines. 

A  summarization  ot  SAR  derivations  and  implied  results  is  given  here.  First, 
isodoppler  lines  and  isorange  lines  are  drawn  tor  a  tlat  surtace.  It  the 
surtace  includes  hills,  isodoppler  lines  move  toward  the  Y  axis  and  isorange 
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lines  move  toward  the  origin.  A  depression  in  the  surtace  causes  isodoppler 
lines  to  move  toward  the  Y  axis  and  isorange  lines  to  move  trom  the  origin. 
Conventional  radar  range  changes  as  a  tunction  of  terrain  elevation,  although 
azimuth  mapping  is  not  attected.  In  synthetic  aperture  radar,  range  and 
azimuth  mapping  is  a  tunction  ot  terrain  elevation.  SAR  mapping  is  discussed 
in  later  text.  Another  important  tunction  that  requires  emphasis  is  nonlinear 
azimuth  mapping  in  SAR.  The  isodoppler  lines  are  hyperbolas. 

It  SAR  processing  is  pertormed  without  correction  tor  nonlinear  tactors,  a 
nonlinear  sweep  must  be  used  tor  the  SAR  display.  An  example  ot  nonlinear 
mapping  is  shown  in  Figure  3.1-4,  whete  the  shaded  area  is  bounded  by  two 
isorange  lines  and  by  two  Isodoppler  lines.  One  isodoppler  line  azimuth 
varies  trom  0  to  19  degrees  and  the  other  isodoppler  line  varies  trom  23  to  28 
degrees. 

3.5  FOCUSED  ARRAY 

The  theory  ot  tocused  array  is  included  as  a  means  ot  understanding  simulation 
considerations.  In  early  SAR  systems,  the  angle  ot  view  was  untocused  and  90 
degrees  trom  the  velocity  vector.  An  untocused  SAR  is  equivalent  to  a  pinhole 
camera  in  optics.  Figure  3.5-1  is  a  diagram  ot  an  untocused  array.  Point  P 
is  the  reterence  point  (the  azimuth  center  ot  a  picture),  line  C  is  the  path 
length  at  the  start  ot  data  collection,  and  line  A  is  the  path  length  at  the 
end  ot  data  collection.  Line  B  is  the  entire  length  ot  data  collection  (the 
aperture  width).  The  law  ot  cosines  is  used  to  derive  the  equation  tor  an 
untocused  array.  This  relationship  is  given  in  equation  3.5-1. 

Af*=Djl+  CS~  2BB  Ccxot4>  (3.5.1) 

Where 

A.  =■  R-XK 
A  p 

B  =  data  collection  distance 

B 

K 

P 


=  phase  shitt  limit 
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Angle  4>  is  the  squint  angle. 

Let  A  *  ft  -X»K  and  C  =  R  in  Equation  3.5.1  to  give: 

A  p  C 

R*~  2  I?  X  Kp  +  (X  K,)2  **Bi  +  Rs-  2Bb  R  cos(^»)  (3.5.2) 


The  solution  ot  (the  untocused  path  distance)  is  given: 
Bb  =  ««*(♦)  +  -  2  \  KP/R  +  (X ■Kp/R)*]'2 


(3.5.3) 


For  veritication,  let  the  squint  angle  (♦)  be  90  degrees  and  delete  the  term 
Equation  3.5.3  now  becomes:  (X*K/R)/2. 

Bb  =  \/2\KpR  (3.5.4) 


Kp  represents  the  amount  ot  phase  shift  that  can  be  tolerated.  Forty-five 
degrees  phase  shift  is  represented  by  a  Kp  ot  one-eighth.  At  90  degrees  the 
untocused  distance  is  doubled  because  of  the  symmetry  around  90  degrees. 
Putting  these  values  in  Equation  3.5.4  gives  the  equation  ot  an  untocused 
array  at  90  degrees: 

Bb  =x/xlT  (3.5.5) 


This  equation  is  given  in  "Introduction  to  Synthetic  Array  and  Imaging  Radars" 
by  S.A.  Hovanessin  (page  23). 


n< 
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By  establishing  a  certain  confidence  in  Equation  3.3.3,  we  turn  to  the 

implications  when  the  squint  angle  is  not  90  degrees.  A  realistic  value  tor 

2*X*K  /R  is  0.0000004.  After  the  cosine  of  4>  exceeds  this  value,  the 
P 

unfocused  distance  approaches  zero.  SAR  will  not  function  in  an  unfocused 
mode  when  the  squint  angle  is  not  90  degrees. 


mi 


Let  angle  4>  be  zero  degrees  and  solve  tor  B  in  Equation  3.5.3.  Travel 

B 

distance  tor  an  unfocused  array  is  X«K  .  This  distance  gives  a  SAR 

P 

aperture  ot  approximately  one-eighth  of  an  inch.  An  unfocused  array  varies 

1/2 

from  zero  at  zero  degrees  to  (X*Kp»R)  at  90  degrees. 

An  important  question  tor  the  simulation  design  engineer  is  a  satisfactory 
solution  it  focusing  tails. 


3.6  SAR  AZIMUTH  RESOLUTION 

There  are  actually  three  resolutions  to  consider  in  simulating  SAR  azimuth 
resolution.  These  resolutions  are:  11  aircraft  display  pixel  resolution, 
2)  aircraft  SAR  resolution,  and  3)  data  base  resolution.  To  make  use  ot  the 
SAR  resolution,  the  pixel  resolution  should  be  at  least  twice  that  ot  the  SAR 
resolution.  This  relationship  is  from  the  Nyquist  Sampling  Theorem.  Since 
simulation  is  to  be  performed  digitally  an  equation  is  needed  to  relate  the 
three  resolutions.  This  equation  will  not  be  an  exact  equation  but  a 
preferred  or  "best  tit"  solution.  Pixel  resolution  will  be  the  pixel 


resolution  ot  the  aircralt.  The  SAR  resolution  will  then  be  modified  to  be 
compatible  with  the  pixel  resolution.  The  data  base  resolution  will  be  made 
to  tit  the  above  resolutions.  The  equation  relating  the  three  resolutions: 


f^DS  _  Nl  RsAK  _  A'2  Rdb 

Kk  Kk  Kk 


Ni>\  N2>2 


(3.6.11 


--VV-1 

'.->1 
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Where 


=  aircratt  display  resolution 
=  synthetic  aperture  radar  resolution 
=  data  base  resolution 


i'- 

•V  ^■---1 

*  c;>:- 


=  selectable  constant  tor  SAR  resolution 
N  =  integers 


The  tactor  K  is  directly  related  to  the  computing  required.  The  ideal 
K 

value  ot  K  would  be  equal  to  one  to  keep  the  computing  load  down,  but  it 
K 

should  be  higher  it  the  response  distribution  ot  a  point  is  to  be  simulated 

with  greater  detail  when  N1  =  2.  The  data  base  should  be  equal  to  the  SAR 

resolution  tor  high  quality  simulation.  In  the  tirst  example,  resolution  is 

considered  to  be  the  size  ot  one  pixel,  but  aircratt  SAR  resolution  and  data 

base  resolution  are  considered  as  two  pixels.  The  minimum  value  ot  K  is 

K 

one  and  Ml  and  N2  have  minimum  values  ot  two.  Assume  that  the  antenna  6  dB 
beamwidth  is  two  pixels,  the  antenna  null  to  null  is  tour  pixels,  and  the 
tirst  side  lobes  are  down  20  dB  and  spaced  three  pixels  trom  the  center  ot  the 
main  lobe.  Figure  3.6-1  is  a  plot  ot  the  response  ot  a  corner  ret  lector 
(point).  Figure  3.6-2  shows  two  corner  ret  lectors  spaced  at  a  distance  ot 
eight  pixels.  Figures  3.6-3,  3.6-4,  3.6-5,  and  3.6-6  illustrate  what  happens 
as  the  two  corner  ret  lectors  (points)  move  together  to  the  limit  ot 
resolution.  The  numbers  chosen  tor  this  example  are  the  minimum  numbers  that 
can  be  used  to  simulate  the  amplitude  distribution  ot  a  point.  The  above 
example  maps  the  point  response  into  live  locations  in  azimuth,,  at  the  same 
slant  range. 

Resolution  in  a  SAR  system  can  be  greater  than  two  pixels.  The  two-pixel 
example  is  the  ideal.  Figure  3.6-7  is  an  example  ot  the  amplitude 
distribution  ot  a  point  when  the  SAR  resolution  is  six  pixels.  The  point  now 
maps  into  23  pixel  locations.  This  should  not  happen  in  actuality  as  the 
optimum  SAR  resolution  in  an  aircratt  is  two  pixels.  Reter  to  Figures  3.6-8 
and  3.6-9  tor  examples. 
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■SOLUTION 


"To  quantify  the  HRR  resolution  in  range  and  azimuth"  was  an  important  objective  for  the  evaluation  flight  test.  An  array  of 
nine  corner  reflectors  was  erected  on  Edwards  AFB  range  PB-4  to  support  these  measurements.  The  final  configuration  of  the 
array  is  shown  below. 
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Figure  3.6-9  shows  the  setup  used  in  evaluating  the  F- 15  SAR.  Figure  3.6-9 
is  a  photograph  ot  the  SAR  display  when  viewing  test  targets.  In  the  16-X 
zoom  ot  the  SAR  map.  the  two  targets  that  are  spaced  15  teet  apart  cannot  be 
separated,  but  the  targets  that  are  spaced  20  teet  apart  can  be  identitied  as 
two  separate  targets.  Also  note  that  azimuth  resolution  and  range  resolution 
seem  to  track. 

3.7  RANGE  WALK  AND  ROTATION 


In  the  theoretical  evaluation  ot  SAR,  it  was  assumed  that  data  would  be 
collected  along  an  arc  and  processed.  Upon  closer  examination  ot  the  details 
ot  SAR  tunctioning,  SAR  processing  becomes  more  complex.  Two  processes,  range 
walk  and  rotation,  tend  to  make  SAR  processing  more  complex.  Figures  3.7-1 
and  3.7-2  are  diagrams  ot  the  two  processes  and  trom  Figure  3.7-2,  equations 
tor  range  walk  and  rotation  can  be  developed.  The  equation  tor  range  walk  is 
the  trigonometric  cosine  law. 


Rw  =*R i  -  s/R?  +Rc-  2  Rx  Pc  cot*) 


13.7.1) 


Where 

R  =  range  walk 

w 

R^  =  range  at  position  one 

The  equation  tor  rotation  is  the  trigonometric  cosine  law  and  the  sine  law. 

(3.7.2) 

(3.7.3) 

Where 

R^  -  range  at  position  two 

-  squint  angle  at  position  two 
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For  better  understanding,  an  example  ot  range  walk  and  rotation  is  indicated 
in  the  tollowing  parameters. 

Range  10  nautical  miles 

Wave  length  3  centimeters 

Squint  angle  45  degrees 

Resolution  1  cycle  in  20  teet  (10-toot  pixel) 

Number  ot  pixels  512 

Path  length  can  be  solved  by  using  Equation  3.2.23  (Section  3.2).  It  the  path 
length  is  0.052180  nautical  miles,  the  range  walk  is  111.41  teet,  and  this 
distance  is  11  pixels.  It  is  apparent  that  range  walk  makes  S/^R  processing 
much  more  ditticult.  Designers  ot  airborne  SAR  systems  have  solved  range  walk 
in  a  very  elegant  manner.  The  pulse  repetition  rate  (PRF )  is  programmed  to 
remove  range  walk. 

Rotation  problems  still  exist,  however,  and  an  example  ot  rotation  is  a 
continuation  ot  the  range  walk  problem.  The  solution  tor  squint  angle/two  is 
45.2107  degrees.  At  a  range  ot  10  nautical  miles,  512  ten-toot  pixels  cover 

an  arc  ot  4.825  degrees.  When  these  numbers  are  put  into  equations,  it  is 

tound  that,  at  their  ends,  the  arcs  are  displaced  trom  each  other  by  1.2 
pixels.  The  designer  ot  the  aircratt  SAR  can  handle  this  problem  in  one  ot 
two  ways.  The  loss  ot  resolution  at  the  edge  ot  the  picture  can  be  accepted 
or  the  problem  can  be  resolved.  In  a  spotlight  SAR  system,  this  migration  is 
corrected  so  that  very  high  resolution  can  be  obtained. 

In  the  design  ot  a  SAR  simulator,  there  must  be  a  clear  determination  ot  what 

is  happening  in  the  actual  aircralt  SAR  system.  It  the  problem  has  been 

eliminated  in  the  aircralt  SAR,  there  is  no  longer  a  problem  in  the  simulator 
SAR.  When  the  rotation  problem  does  exist  in  the  aircratt,  there  is  a  simple 
solution  available  tor  simulator  design:  a  small  amount  ot  weighted  averaging 
at  the  edges  ot  the  picture  will  give  the  right  amount  ot  detocusing. 
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There  is  an  additional  problem  caused  by  rotation.  In  some  SAR  modes,  the 
display  is  made  up  ot  two  or  more  sector  scans.  Figure  3.7-3  shows  how  two 
sector  scans  are  brought  together  to  form  a  single  picture.  Since  the  two 
scans  are  generated  trom  two  locations,  the  two  sectors  will  not  have  a 
pertect  tit.  In  the  aircratt  SAR,  the  problem  can  either  be  ignored  or  there 
can  be  a  small  amount  ot  overscan  and  scan  merger.  Merging  is  accomplished 
using  a  weighted  average  ot  the  pixels  that  overlap.  The  simulator  design 
should  replicate  the  aircratt  SAR  as  a  solution  to  this  problem. 

3.8  PULSE  COMPRESSION 

As  stated  previously,  radar  has  two  dittering  resolutions  tor  range  direction 
and  azimuth  direction.  The  range  resolution  ot  conventional  pulse  radar  is 
given  by  the  tol lowing  equation: 

AR  «  c  t/2  (3.8.1) 


Where 

C  =  velocity  ot  light  v 

R  =  range  resolution  in  meters 

Figure  3.8-1  presents  a  graphic  analysis  ot  the  derivation  ot  the  above 
equat ion. 

Radar  return  signal  power  varies  as  the  inverse  tourth  power  ot  range.  Radar 
power  is  a  function  ot  the  pulse  area  and  radar  range  resolution  varies 
inversely  with  pulse  width.  A  problem  occurs  it  a  high  powered  radar  with 
high  range  resolution  is  desired.  In  this  circumstance,  the  best  solution 
would  be  to  build  a  radar  with  zero  pulse  width  and  infinite  peak  power.  With 
conventional  pulse  radar,  some  type  ot  compromise  must  be  made  between  range 
and  resolution. 
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Figure  3.8-1  RADAR  RESOLUTION  DUE  TO  PULSE  WIDTH 


N AVTRAEQU 1PCENT  83-C-0068 


Pulse  compression  allows  pulse  length  to  be  independently  varied  trom  range 
resolution.  Using  this  technique,  range  resolution  becomes  a  tunction  ot  the 
bandwidth  ot  the  radar  transmitter  and  receiver.  * 


Where 

B  =  bandwidth 

W 

This  equation  implies  that  peak  pulse  amplitude  can  be  kept  to  an  optimum 
level  and  range  resolution  can  be  increased,  it  a  way  can  be  bound  to  increase 
the  bandwidth  without  increasing  the  peak  power.  Fortunately  there  are  many 
ways  to  accomplish  this.  First  attempts  showed  a  modest  gain  in  range 
resolution,  with  advances  in  hardware  capabilities,  significant  increase  in 
range  resolution  tol lowed.  In  order  to  understand  what  is  to  be  done  to 
simulate  radars  ot  higher  resolution,  two  examples  ot  techniques  tor  increased 
range  resolution  are  included. 

The  tirst  example  will  be  linear  frequency  modulation  radar,  commonly  referred 
to  as  chirp  radai .  The  chirp  radar  principle  consists  ot  linearly  modulating 
the  carrier  frequency  over  the  length  ot  the  transmitted  pulse  and  putting  a 
matched  filter  in  the  receiver  that  will  separate  returns  in  time,  as  a 
tunction  ot  carrier  frequency.  Figure  3.8-2  illustrates  the  concept  ot  chirp 
radar.  The  trequency  at  the  beginning  ot  the  pulse  is  delayed  until  the  end 
ot  the  pulse;  all  other  frequencies  are  proportionally  delayed.  This  delay 
causes  the  return  pulse  to  Increase  in  height  and  decrease  in  width.  The  cost 
tor  achieving  this  increase  in  range  resolution  is  a  more  complex  radar  and  a 
larger  bandwidth.  There  are  additional  types  ot  carrier  modulation;  however, 
tor  each  type,  a  filter  must  be  designed  to  match  the  modulation.  Chirp 
modulation  has  an  advantage  over  other  methods  because  any  mismatch  between 
the  modulator  and  t liter  causes  only  a  constant  range  error.  One  disadvantage 
ot  basic  chirp  radar  is  the  generation  of  side  lobes  in  the  range  direction. 
These  side  lobes  have  the  same  tunction  as  a  lens  (sinx/x)  and  can  be 
controlled  in  the  same  way.  Uncontrolled  side  lobes  are  down  only  13.2  dB. 
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It  a  chirp  radar  system  is  to  be  simulated,  the  aperture  response  ot  the  chirp 
wavetorm  must  be  determined  and  a  beamspread  algorithm  can  then  be  used  to 
simulate  the  aperture  response. 

A  second  method  divides  the  pulse  into  a  number  ot  bins  and  the  phase  or 
frequency  is  shitted  trom  bin  to  bin.  Although  many  ways  ot  achieving  this 
have  been  proposed  and  tested,  this  report  only  introduces  the  phase  shitted 
Barker  coded  technique.  This  method  seems  to  be  most  favored  by  modern  radar 
designers.  Barker  codes  were  developed  tor  the  communication  industry;  Barker 
codes  tor  2,  3,  4,  5,  7,  11,  and  13  elements  have  been  discovered  and  searches 
tor  larger  codes  continue.  Barker  codes  can  be  combined,  or  combined  with 
other  codes,  to  get  higher  resolution  and  lower  side  lobes.  The  Hughes  F-15 
demonstration  radar  has  a  programmable  pulse  compression  that  has  a  limit  ot 
169  (13  times  13).  The  code  in  the  P-15  demonstration  radar  could  be  a  nested 
Barker  code.  Pulse  compression  codes  tor  radar  on  the  F/A-18  aircraft  are  not 
known  at  this  time.  Two  references  tor  Barker  codes  tor  radar  use  are 
"Principles  ot  High-Resolution  Radar"  by  August  W.  Rihaczek,  and  "Radar 
Handbook"  by  Merrill  I.  Skolnik. 

Figure  3.8-3  is  a  diagram  ot  autocorrelation  ot  a  single  pulse  and  a 
13-element  Barker  coded  pulse.  By  phase  coding  13  pulses,  resolution  is 
increased  by  13  and  12  side  lobes  with  amplitude  ot  1/13  are  generated.  To 
simulate  a  Barker  code  ot  13  requires  a  beamspread  algorithm  that  operates  in 
the  range  direction.  The  designer  ot  high  resolution  radar  systems  generally 
matches  range  resolution  with  azimuth  resolution.  It  the  radar  being 
simulated  does  not  have  equal  range  resolution  and  azimuth  resolution,  the 
simulator  designer  must  take  this  into  account.  Additionally,  it  range  and 
azimuth  side  lobes  can  be  seen  in  the  actual  SAR  system,  the  side  lobes  must 
be  simulated. 

There  is  one  area  ot  pulse  compression  simulation  that  requires  particular 
attention.  Solution  is  needed  tor  pulse  compression  simulation  when  the  radar 
overloads  or  is  jammed.  Figure  3.8-4  is  a  SAR  photograph  that  is  only 
overloaded  by  6dB.  The  two  photographs  indicate  that  when  SAR  overloads,  the 
results  are  quite  dramatic.  A  cross  is  formed  around  the  target;  the  azimuth 


It  resolution 
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portion  consists  ot  a  smear  that  is  one- third  ot  the  picture.  The  range 
portion  consists  ot  a  smear  ot  one-tourth  ot  the  picture,  plus  dots  that  cover 
the  full  picture.  The  photograph  was  made  trom  subpictures  and  only  the 
center  subpicture  was  overloaded  in  the  azimuth  direction.  Since  the  range 
saturation  ettect  extended  across  the  picture,  the  radar  pulse  covered  more 
than  one-halt  the  picture. 

The  saturation  ettect  in  the  range  direction  will  cover  twice  the  pulse  width 
distance.  Figure  3.8-4  illustrates  why  overloads  and  jamming  must  be 
simulated.  Simulation  algorithms  tor  overloads  and  jamming  must  be  custom 
titted  to  each  combat  radar. 
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4.0  RADAR  SIMULATION 

4.1  RADAR  EQUATION  MECHANIZATION 

The  basic  radar  equation  was  formulated  in  the  early  history  of  radar.  Many 
modifications  and  retinements  to  the  basic  equation  have  been  made,  but  the 
basic  equation  remains  the  best  starting  point  for  examining  radar 
simulation.  A  good  representative  radar  equation  is  the  one  developed  by 
D.E.  Kerr  and  "The  Radar  Handbook"  by  Merrill  Skolnlk  also  has  an  excellent 
discussion  of  radar  equations.  The  radar  equation  given  by  Kerr  is  included 
here: 

Pr  Gta'Gra'Tcs^PftPfr 

Pt  "  «•*)•■*«  l4,1’l) 


Adding  the  effects  ot  atmospheric/weather  attenuation  and  radar  receiver  gain 
functions  to  (Gc)  Equation  4.1.1  generates  Equation  4.1.2.  The  terms  in 
Equation  4.1.2  are  also  regrouped. 


Pr 

Pt 


Gta'Gra^Gc 
(4t r)*R4 


(4.1.2) 


Where 

=  transmitted  power 

Pr  =  received  power 

G  =  transmitter  antenna  gain 

TA 

G  *  receiver  antenna  gain 

RA 

P  =  pattern  factor  along  the  transmitting  path 

FT 

PpR  =  pattern  tactor  along  the  receiving  path 

T  =  target  cross  section 

Gc  =  gain  function  due  to  gain  control,  atmospheric/weather  attenuation 

X  =  wavelength 

R  =  range 
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Equation  4.1.2  now  consists  ot  three  basic  parts.  Part  one:  the  gain  tactor 
due  to  range  attenuation,  receiver  gain,  atmospheric/weather  attenuation,  and 
receiver  gain  control.  Part  two:  the  attenuation  due  to  the  pattern  ot  the 
radar  antenna  (this  tactor  varies  Irom  one  to  zero).  Part  three:  the 
retlectivity  ot  the  target. 

In  part  one  ot  Equation  4.1.2,  it  the  gain  control  ot  the  radar  receiver  were 
pertect,  this  part  ot  the  equation  will  be  unity  until  the  receiver  gain  limit 
is  reached  (at  this  time  we  will  not  cover  the  signal- to- noise  problem).  It 
is  obvious  that  the  gain  tunction  (Gc)  varies  with  radar  systems.  Since 
this  study  is  not  directed  to  a  single  radar  system,  the  gain  tunction  ot  part 
one  will  be  considered  to  be  unity  multiplied  by  an  unknown  radar  system 
tunction. 

Part  two  ot  Equation  4.1.2  modities  the  returns  tor  target  radars  and 
synthetic  aperture  radars  and  is  also  responsible  tor  the  beamspread  in  ground 
mapping  real  aperture  radar.  This  tunction  modities  the  return  relative  to 
its  position  in  the  beam.  Since  ground  mapping  real  aperture  radar  returns 
are  the  sum  ot  all  targets  within  the  beam,  the  tunction  spreads  out  the 
returns.  In  tire  control  radar,  the  position  ot  the  target  in  the  antenna 
pattern  determines  the  return  level.  SAR  has  two  apertures  that  determine 
performance.  The  first  aperture  is  the  real  antenna  aperture  which  determines 
the  tootprint  (the  pattern  tactor)  ot  the  return.  The  second  aperture  (the 
SAR  aperture)  determines  the  resolution.  The  upper  lett  section  ot  Figure 
4.1-1  is  an  example  ot  a  SAR  picture  and  the  ettects  ot  an  antenna  pattern. 
Note  the  radar  antenna  was  positioned  11  times  to  obtain  the  SAR  picture.  The 
scalloped  window  blind  ettect  is  due  to  antenna  pattern. 

Part  three  ot  Equation  4.1.2  (the  radar  target  cross  section)  is  a  very 
pertinent  tunction  tor  radar  simulation.  The  radar  target  cross  section 
consists  ot  six  variables.  Hence  the  radar  target  cross  section  can  be 
expressed  as: 


4-2 
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TCS — H 'W -Agp-O-T/fgp-Sg  (4.1.3) 

Where 

H  =  height 

W  =  width 

A  =  aspect 

SP 

O  =  orientation 

S  =  shadow  effect 

E 

T  *  target  reflectivity 

REF 

The  functions  that  make  up  the  radar  target  cross  section  are  analyzed  in 
later  text.  First,  an  evaluation  is  made  of  how  the  components  of  the 
equation  determine  the  range  target  cross  section.  Height  is  modified  by  the 
shadow  effect.  For  example,  it  one  building  blocks  90  percent  of  another 
building,  the  shadow  effect  would  be  0.1.  Thus  the  effective  height  would 
only  be  ten  percent.  It  a  target  were  wider  than  a  pixel,  then  the  width 
would  be  the  width  of  a  pixel.  It  the  target  is  less  than  the  width  of  a 
pixel,  the  width  is  then  less  than  a  pixel. 

Aspect  is  the  angle  at  which  the  radar  beam  intercepts  the  target  or  the 
grazing  angle.  The  aspect  angle  is  a  function  of  the  depression  angle  and  the 
face  angle  of  the  target,  with  respect  to  the  horizon. 

Orientation  is  a  function  of  the  azimuth  angle  that  the  radar  beam  intercepts 

in  the  target.  An  example  of  this  function  could  be  mapping  that  includes  an 

outdoor  movie  screen.  Because  the  front  of  the  screen  slants  downward,  the 

reflected  light  is  directed  toward  the  earth  and  little  energy  is  reflected. 

At  a  90  degree  angle,  the  cross  section  is  small  and  little  energy  is 

reflected.  The  back  side  of  the  screen  is  tilted  toward  the  aircraft  and  is  a 

textured  surface;  thus,  the  radar  return  is  strong.  T  is  the  actual 

REF 

reflectance  of  the  target.  Every  surface  material  has  a  different  reflectance 
level;  tor  instance,  metal  is  highly  reflective  and  grass  is  almost  non- 
retlective.  In  simulated  radar,  the  reflectance  value  is  derived  from  the 
data  base  and  is  not  calculated. 
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4.2  ATTENUATION  RESULTING  FROM  ORIENTATION 

DMA  Level  I,  Level  II,  Level  V,  and  Level  X  data  bases  have  provisions  tor 
coding  radar  returns  as  a  tunction  of  azimuth.  Point  features  have  the 
maximum  return  specitied  on  11.25  degree  increments  trom  true  north,  or 
spec  it  led  as  omnidirectional.  Linear  teature  coding  ot  returns  can  be 
unidirectional,  bidirectional,  and  omnidirectional.  Unidirectional  is  detined 
as  a  retlective  surtace  on  the  right  side  ot  a  line  trom  Point  1  to  Point  2. 
Bidirectional  is  detined  as  a  retlective  surtace  on  both  sides  ot  the  line 
trom  Point  1  to  Point  2.  Omnidirectional  is  detined  as  a  surtace  with  equal 
reflectivity  on  all  sides. 

Simulation  ot  orientation  is  calculated  to  the  nearest  11.25  degree 
increment.  The  teature  identification  code  number  (FICN)  and  orientation  are 
then  used  to  obtain  a  reflectance  multiplier  trom  a  look-up  table.  A  large 
number  ot  teatures  will  have  identical  orientation  tunctions. 

4.3  ASPECT 

One  ot  the  six  tunctions  that  determine  the  cross  section  ot  the  radar  is  the 
aspect  or  grazing  angle.  A  simulation  algorithm  can  be  derived  that  is  accu¬ 
rate  and  easily  implemented.  The  grazing  angle  return  tunction  is  determined 
by  the  grazing  angle  (or  the  angle  ot  incidence),  radar  wavelength,  radar  beam 
polarization,  amount  ot  moisture  on  the  target  surtace,  and  by  surface  vegeta¬ 
tion  and  cultural  teatures.  A  simplified  algorithm  tor  aspect  is  a  determina¬ 
tion  ot  grazing  angle  and  target  type  and  a  look-up  table  entry  of  the  proper 
value  ot  the  multiplier  in  order  to  obtain  the  reflectance  value  ot  the  target. 

An  lntinite  number  ot  look-up  tables  could  be  generated  tor  each  set  of  aspect 
conditions,  but  only  a  tew  tables  are  actually  needed.  Simulated  seasonal 
ettects  are  generated  trom  aspect  tables.  For  instance,  a  field  ot  summer 
grass  retlects  some  return  at  all  grazing  angles.  In  winter,  it  heavy  snow 
had  covered  the  same  area,  forward  scatter  of  the  return  would  be  increased 
and  back  scatter  return  would  be  decreased.  In  some  circumstances  the  area 
could  disappear  trom  the  radar  screen. 


4-5 


NAVTRAEQUIPCENT  83-00068 


Figure  4.3-1  is  a  generalized  curve  versus  grazing  angle.  The  generalized 
curve  changes  very  little  in  the  center  medium  grazing  angle.  At  a  high 
grazing  angle  the  curve  has  a  high  return,  reterred  to  as  clilt  ettect  or 
vertical  incidence.  As  the  grazing  angle  approaches  zero,  return  approaches 
zero.  This  grazing  angle  ettect  is  responsible  tor  hilltop  and  distance 
lading  and  tor  loss  ot  resolution  at  long  ranges.  In  simulating  radar,  all 
grazing  angles  are  ot  concern,  but  high  and  low  grazing  angles  should  receive 
greater  emphasis  and  attention. 

Terrain  aspect  effects  are  generated  from  the  elevation  data  base.  Any  data 
base  limitations  restrict  aspect  effects.  For  example;  consider  simulation  ot 
a  flat  plane,  a  300- toot  vertical  clitt  and  continuing  flat  plain,  viewed  at  a 
low  angle.  The  Defense  Mapping  Agency  Aerospace  Center  (DMAAC)  elevation  data 
base  is  on  300-toot  posts.  A  300- toot  cliff  will  give  a  45  degree  aspect 
angle  with  this  post  spacing.  If  the  cultural  data  base  has  a  cliff  feature, 
the  clitt  is  treated  properly;  otherwise,  the  return  is  treated  as  a  45  degree 
aspect  angle.  Most  man-made  structures  have  horizontal  or  vertical  flat 
surfaces  and  have  the  greatest  returns  when  viewed  at  a  low  grazing  angle 
(perpendicular  or  upright).  Figure  4.3-2  is  a  diagram  ot  a  grazing  angle 
curve  tor  a  teature  with  a  flat  roof,  and  Figure  4.3-3  is  a  diagram  ot  grazing 
angle  curve  tor  a  teature  with  a  peaked  roof.  Conversely,  a  rounded  and 
highly  reflective  root  is  almost  invisible  to  radar. 

An  aspect  algorithm  consists  of  the  sum  ot  the  depression  angle  and  the  slope 
ot  the  terrain.  It  the  target  is  a  cultural  object,  terrain  slope  is  set  to 
zero.  Figure  4.3-4  is  a  vertical  aspect  diagram. 

The  equation  tor  the  slope  in  Figure  4.3-4; 

(4.3.1) 

Where 

o  =  terrain  slope  ot  a  line  connecting  El  &  E2 

R  -  delta  ground  range 

DG 

E2  S>  El  =  two  adjacent  elevation  posts 


£2  -  £1 
a  =tan  — - - 
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Depression  angle  calculation  becomes  more  complex  it  earth  curvature  is 
considered.  Figure  4.3-5  is  a  diagram  ot  geometry  used  to  calculate  the 
depression  angle.  Value  ot  the  earth's  radius  can  be  adjusted  to  accommodate 
the  earth's  atmospheric  retraction;  a  commonly  used  value  is  tour- thirds.  By 
using  simple  geometry,  the  equation  tor  depression  angle  (Equation  4.3.2)  can 
be  obtained. 


V>  =  tan 


,  R  +  REC  -  cos(()  Rec 
sin(0  Rec 


(4.3.2) 


simplitying  equation  4.3.2  gives 


=  tan  1 


WfcsioU) 


(4.3.3) 


Where 

R  =  modi  tied  earth  radius  tor  atmospheric  retraction 
EC 

5  =  ground  range  angle 


The  grazing  angle  tor  terrain  is; 


7  =a  + 


(4.3.4) 


The  grazing  angle  tor  culture  objects  is: 


7  —i> 


(4.3.5) 


Where 

Y  =  grazing  angle 


The  SARS  used  in  combat  aircratt  have  the  capability  ot  generating  images  at 
large  depression  angles.  Most  SAR  systems  are  limited  to  65  degrees  or  less. 
This  means  that  a  look-up  table  trom  0  to  65  degrees  tor  culture  targets  would 
suttice.  A  terrain  look-up  table  trom  0  to  65  degrees  would  match  the  culture 
look-up  table.  Computer-generated  images  have  shown  that  one-quarter  degree 
increments  are  sutticient  tor  the  look-up  table. 


,  & 


V  W  V  V  Y»  T>  Y 


'*  o*  *Y  *  *  * 
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Equation  4.3.3  can  be  simplified  at  short  ground  ranges  by  ignoring  the 
curvature  of  the  earth.  Equation  4.3.6  is  this  simplification. 


V>  “tan' 


(4.3.6) 


It  has  been  brought  out  that  reflectance  is  a  function  of  the  grazing  angle 
and  the  orientation  of  the  target.  The  grazing  angle  can  be  calculated  tor 
all  targets  but  orientation  can  only  be  calculated  when  the  orientation  is  in 
the  data  base.  Level  I  and  Level  II  data  bases  code  the  orientation  for  point 
and  linear  features.  Level  V  data  bases  include  orientation  for  point, 
linear,  and  aerial  features.  These  features  can  be  randomly  selected  to  have 
either  north/south  orientation  or  east/west  orientation.  This  function  is 
called  the  Homestead  Law  Effect. 

Farms  in  the  mid-West  and  also  the  central  valley  of  California  were  laid  out 
on  north/south,  east/west  grids.  Farmers  plow  parallel  to  these  grid  lines 
and  radar  returns  are  stronger  at  the  cardinal  points  of  the  compass.  To 
determine  the  number  of  aspect  and  orientation  look-up  tables  that  are 
required,  each  feature  specified  in  the  DMAAC  data  base  must  be  analyzed. 

4.4  SLANT  RANGE  AND  DEPRESSION  EFFECTS 

Two  important  parameters  in  radar  simulation  are  slant  range  and  depression 
angle.  Slant  range  is  used  to  map  pixels  from  the  ground  range  scan  into  the 
display  scan.  The  depression  angle  is  used  to  calculate  both  the  aspect 
effect  and  the  shadow  effect.  Figure  4.4-1  is  a  diagram  used  to  derive 
equations  tor  slant  range  and  depression  angle. 

The  earth's  atmosphere  causes  a  certain  amount  of  refraction  of  the  radar 
signals  and  a  simple  method  of  simulating  this  refraction  is  to  increase  the 
earth  radius  during  simulation  calculations.  A  value  most  often  used  is 
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tour-thirds.  In  the  equations,  the  modification  ot  the  earth's  radius  will  be 

referred  to  as  R  .  The  range  angle  is  related  to  range  by  Equation  4.4.1. 
EC 

Figure  4.4-2  is  a  diagram  tor  deriving  equations  in  slant  range/top  of 
structure  mapping. 


Ra 


Where 

K  =  earth  radius  modification  factor 
E 


(4.4.1) 


The  law  of  cosines  can  be  used  to  solve  for  slant  range  in  Figure  4.4-1. 

RsL  —  J(^£C  +  -A)  +  (^£C  +  -  {ReC  +  +  £)c°s(0]  ( 


(4.4.2) 


By  using  the  law  ot  sines  and  the  slant  range,  the  depression  angle  can  be 
derived. 


-i  -  ,  (gec  + 

V-c  =6in  8in(0 - 5 — 

ksl 


(4.4.3) 


Therefore 


=cos_1  pin(() 


(Rec  +  E) 


Noting  that  ^  +  T  -  90‘ 


(4.4.4) 


Where 


=  depression  angle  complement 


By  letting  the  line  E  remain  parallel  to  line  A,  more  simple  equations  can  be 
derived,  although  with  a  slight  error.  The  depression  angle  is  given  by 
Equation  4.4.5  and  the  slant  range  is  given  by  Equation  4.4.6. 


tj>  =*  tan” 1 


L  -  E  + 


Rec\  1  -  cos 


—  «n(0l 


(4.4.5) 
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WHERE 

R<..T  *  SLANT  RANGE  AT  THE  TOP  OF  THE 
STRUCTURE 

R-,  =  SLANT  RANGE  AT  THE  BOTTOM  OF 

THE  STRUCTURE 
UJ  =  DEPRESSION  ANGLE 


Figure  4.4-2  SLANT  RANGE  MAPPING  DIAGRAM:  TOP  OF  STRUCTURE 
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Where 

E  =  elevation 


Rsl 


=*Rec 


.  s»n( Q 
sin  (V) 


(4.4.6) 


The  top  of  a  building  will  map  toward  the  aircraft.  An  equation  will  be 
derived  to  show  what  happens  when  mapping  high  structures.  The  equation  will 
express  the  error  from  the  bottom  to  the  top  of  the  building.  Figure  4.4-2  is 
the  diagram  that  will  be  used. 

One  task  is  to  find  the  difference  in  lengths  between  lines  R  and  R 

SL  SLT 

By  the  law  ot  cosines.  Equation  4.4.7  can  be  derived: 

Rslt  —y/E*  +  Rsl  -  2  E  RSL  cos[\l>e)  (4.4.7) 


Expanding  Equation  4.4.7  by  using  the  binomial  series  gives: 


Rslt  —Rslt 


1  + 


E2 
2  Rsl 


2EEsx,cob(V'c) 

2 rIl 


(4.4.8) 


2  2 

Since  E  /RSL  is  much  less  than  one  and  cos(Tc)  is  equal  to  sin(T) 
Equation  4.4.8  can  be  written: 

Rsl  «  Rsl  -  (4.4.9) 


Taking  the  ditterence  between  the  slant  range  at  the  bottom  ot  a  structure  and 
the  top  ot  a  structure  will  give  the  slant  range  mapping  error  at  the  top  ot  * 


structure. 
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Rfu  ^  -  Esin(^) 


(4.4.10) 


Where 

R  =  range  torward  mapping  due  to  elevation 
FM 

Equation  4.4.10  shows  that  the  error  in  slant  range  mapping  at  the  top  ot  a 
structure  is  directly  related  to  the  height  ot  the  structure  and  the 
depression  angle. 

To  get  an  idea  ot  this  mapping  ot  tall  structures,  two  examples  will  be  given. 
Example  One 


Depression  Angle 
Range 

Structure  Height 
Pixel  Size 
Resolution 

Slant  Range  Error 
Pixels  Mapped  into 

Example  Two 

Depression  Angle 
Range 

Structure  Height 
Pixel  Size 
Resolution 

Slant  Range  Error 
Pixels  Mapped  Into 


15  Degrees 
10  Nautical  Miles 
1200  Feet 
10  Feet 
20  Feet 

310  Feet 
31  Pixels 


15  Degrees 
10  Nautical  Miles 
150  Feet 
50  Feet 
100  Feet 

38.8  Feet 
2  Pixels 
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4.5  SHADOWS 

The  radar  shadow  algorithm  is  one  of  the  most  important  algorithms  tor  SAR 
simulation.  The  depression  angle  of  SAR  is  quite  low;  thus,  in  any  SAR  radar 
display  there  will  be  large  shadows.  For  proper  simulation,  the  shadow 
algorithm  must  realistically  duplicate  the  actual  radar.  Three  shadow 
algorithms  will  be  presented  and  defined.  Choice  of  algorithm  is  based  on 
complexity  versus  realism. 

The  first  algorithm  that  will  be  presented  is  the  on/off  shadow  algorithm. 
This  was  the  first  and  most  frequently  used  algorithm  in  radar  simulation. 
Figure  4.5-1  is  a  diagram  ot  this  algorithm.  The  term  on/off  algorithm  is 
used  because  the  target  either  has  a  full  reflectivity  or  has  none.  The 
algorithm  is  simple.  The  algorithm  saves  the  smallest  depression  angle  and 
compares  the  following  return  depression  angle  with  the  saved  smallest 
depression  angle.  when  the  return  depression  angle  is  greater  than  the 
reference  angle,  the  return  reflectance  is  set  to  zero.  When  the  return 
depression  angle  is  less  than  the  reference  angle,  the  return  reflectance  is 
unchanged;  the  return  depression  angle  replaces  the  reference  angle.  At  low 
depression  angles  this  algorithm  is  not  very  realistic. 

The  second  shadow  algorithm  can  be  called  the  partial  shadow  algorithm.  This 
algorithm  modifies  the  radar  return  as  a  function  of  the  degree  of  shadowing 
ot  the  target.  Figure  4.5-2  is  a  diagram  of  this  algorithm.  This  algorithm 
ditters  trom  the  first  shadow  algorithm  by  the  addition  of  a  shadow  function 
that  lies  between  the  on/oft  function.  This  function  modifies  the  target 
return  by  the  amount  of  snadowing.  When  the  target  is  half-shadowed  the 
target  return  is  reduced  by  halt.  An  example  ot  how  this  algorithm  increases 
realism  is  to  look  at  a  village  that  is  observed  at  a  low  grazing  angle.  Only 
the  first  row  of  houses  will  be  fully  reflective.  All  other  houses  will  be 
partially  shadowed  by  the  first  row  of  houses.  The  partial  shadow  algorithm 
greatly  increases  the  realism  ot  radar  simulation. 
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The  third  algorithm  is  used  when  SAR  simulation  is  used  tor  tactical  training 
that  requires  a  high  degree  ot  navigational  accuracy.  Literature  that  is 
available  on  the  use  ot  SAR  in  a  tactical  situation  suggests  that  bridges, 
towers,  and  piers  are  some  ot  the  most  important  radar  features  to  update  the 
inertial  plattorm.  Since  bridges  are  above  terrain,  a  bridge  does  not  shadow 
trom  the  terrain  level  to  the  top  ot  the  bridge,  but  rather  iron  the  bottom  ot 
the  bridge  to  its  top.  This  means  that  tor  proper  shadowing  a  more  complex 
shadowing  algorithm  must  be  used.  Figure  4.5-3  is  a  diagram  ot  this 
algorithm.  The  data  base  must  supply  the  top  and  bottom  bridge  elevations. 
This  algorithm  provides  the  proper  signatures  ot  bridges  and  any  target  that 
has  an  overhang. 

Three  algorithms  tor  the  radar  shadow  have  been  presented.  The  choice  ot 
algorithm  tor  implementation  is  a  matter  ot  cost  versus  realistic  SAR 
simulation.  There  are  more  complex  shadow  algorithms  that  could  have  been 
presented,  but  the  more  complex  algorithms  deal  with  radars  (i.e.,  naval 
warship  signatures)  that  are  not  as  close  to  being  operational  as  SAR  systems. 

One  important  tact  that  should  be  emphasized  is  that  in  all  shadow  algorithms 
tor  SAR  simulation,  shadow  calculations  must  be  performed  from  the  nadir  point 
to  the  end  ot  the  scan.  An  example  ot  this  is  a  one  square  mile  patch  scan  at 
a  40  mile  range.  All  ot  the  SAR  simulation  algorithms  operate  only  on  the  one 
mile  patch  except  the  shadow  algorithm.  Th  shadow  algorithm  operates  on  40 
times  as  much  elevation  data.  Since  tactical  aircraft  are  going  to  fly  as  low 
as  possible,  low  level  SAR  simulation  is  essential. 
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5.0  SAR  ANOMALIES  AND  SPECIAL  EFFECTS 


►v' 

m  , 

i.v 


(:.* 


I? 

r\ 


£ 


5.1  MOVING  TARGETS 

A  surprising  phenomenon  ot  synthetic  aperture  radar  that  does  not  occur  in 
real  aperture  radar  is  the  mapping  characteristics  ot  moving  objects.  In  SAR 
moving  objects  are  mapped  at  a  different  location  than  in  RAR  but  the  moving 
objects  leave  their  shadows.  In  the  simulation  of  SAR,  mapping  of  moving 
objects  is  required.  Moving  targets  of  interest  could  be  trucks,  tanks,  etc. 
This  relocation  ot  moving  objects  in  SAR  is  a  function  of  the  object's 
relative  heading  and  velocity.  The  development  of  an  algorithm  tor  simulation 
is  straighttorward. 

The  Doppler  frequency  is  represented  by  Equation  5.1.1. 


Fds  =2  Vac  F/C 


(5.1.1) 


I 

<v 

f>: 

c 


1  % 


£ 


Let 

Vp  be  the  velocity  along  the  vector  UP  the  aircraft  nadir  point  to  a  point 
on  the  ground. 


The  problem  is  approached  by  calculating  the  Doppler  frequency  shift  tor  a 
moving  object  and  finding  a  stationary  location  on  the  ground  that  will 
generate  the  same  frequency.  Figure  5.1-1  is  a  diagram  of  the  velocity 
vectors  tor  the  moving  object  problem.  To  start,  the  point  ot  reference  will 
be  the  nadir  point.  The  velocity  along  the  direction  ot  the  nadir  point 
to  the  ground  point  P  is: 


Vp  rm  V^c-COsJ©  w) 


(5.1.2) 


Equation  S.1.3  relates  the  moving  object  velocity  to  the  nadir: 
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VUp  ««  V/cr  cos(9r  -  ©*  +  8jt) 


(5.1.3) 


(here 


/  =  Moving  target  velocity  relative  to  the  ground 

/  =  Velocity  between  the  aircratt  nadir  point  and  a  point  on  the  ground 

MP 

extcluding  the  aircratt  velocity 


9  =  Target  heading 

T 


9^=  Aircratt  heading 

9^  =  Angle  between  aircratt  ground  track  and  vector  trom  nadir  to  the 
target 


Equations  5.1.2  and  5.1.3  can  be  combined  to  give  the  total  velocity  with 
respect  to  the  nadir: 


Vi wr  "*  K»Cco8(©>m)  -  ^cr'co* -  ©^  +  ©w) 


(5.1.4) 


Equation  5.1.4  can  now  be  reterenced  to  the  aircratt  by  multiplying  the  cosine 
ot  the  depression  angle: 


VTAC  *cos(V’)  [v^c  co»(®vt)  -  Vcr-co8 {© r  -  ©*  +  ©w)l 


(5.1.5) 
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Where 

V  ^  -  aircratt  velocity  relative  to  a  target  on  the  ground 


The  relative  velocity  between  a  point  on  the  ground  and  the  aircratt  is: 

» 

Vgac  —  Vac  ‘cos(  t)  cos  (8  w )  (5.1.6) 


Where 

VGAC  *  aircratt  velocity  relative  to  to  the  ground 


Equation  5.1.7  is  the  result  ot  equating  Equations  5.1.5  and  5.1.6  into 
Equation  5.1.1. 


VAC ^os(V’)  cos(ew)  —  cosV'^y(c  co8(eVrV)  -  VGT  cos (0 r  -  84  +8yv)j 


(5.1.7) 


8W  —cos  ^os |0 yy |  ~  Vcr/V„ccos(0r-  0*+ 8yv)j 


(5.1.8) 
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Where 

0  3  azimuth  mapping  angle 

H 

Equation  5.1.8  is  a  complete  description  of  azimuth  tor  a  moving  object.  The 
shadow  ot  a  moving  object  is  calculated  using  the  original  azimuth.  The  slant 
range  is  not  influenced  by  the  motion  ot  the  moving  object.  Although  there  is 
some  interaction  between  object  motion  and  some  pulse  compression  algorithms, 
it  is  not  documented  in  this  report.  An  object  moving  from  the  aircraft  maps 
away  from  the  velocity  vector  and  an  object  moving  toward  the  aircraft  maps 
toward  the  velocity  vector.  Moving  objects  in  the  radar  beam  with  their 
ground  location  out  of  the  radar  display  can  map  into  the  radar  display. 
Also,  moving  objects  in  the  radar  display  can  map  out  of  the  radar  display  and 
leave  a  remaining  shadow.  An  example  of  a  moving  object  is  given  below. 


Example:  Moving  Object 


Range 

Aircraft  Heading 
Target  Heading 
Azimuth  Angle 
Aircraft  Velocity 
Target  Velocity 


*  10  nautical  miles 

3  45° 

=  90° 

=  20° 

=  300  knots 
=  10  knots 

=  24.75* 


At  a  range  of  10  nautical  miles,  a  moving  object  would  map  4852  feet  to  the 
right  ot  its  shadow. 
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Equation  5.1.8  is  not  a  complete  answer  when  considering  the  real  world. 
Collecting  data  tor  a  SAR  display  requires  a  finite  amount  ot  time.  A  moving 
object  could  be  traveling  in  a  curved  path.  When  this  occurs  the  signature  is 
many  velocities  and  the  moving  object  is  spread  over  a  number  ot  pixel 
locations  rather  than  one.  A  ditticult  case  tor  simulation  would  be  that  ot  a 
long  train  traveling  on  a  curve,  in  which  the  curvature  is  at  the  midpoint  ot 
the  train. 

Figures  5.1-2  and  5.1-3  are  photographs  ot  a  train  moving  along  a  curved 
track.  For  moving  objects,  the  reflectance  pixel  should  be  spread  over 
several  pixels. 

5.2  NOISE 

The  return  signal  ot  radar  is  attenuated  by  range  to  the  fourth  power.  In 
order  to  simulate  noise  in  a  SAR  picture,  some  type  ot  noise  model  must  be 
developed.  Assume  that  a  target  occupying  the  full  area  ot  a  pixel  is  100 
percent  reflective  and  generates  a  signal  that  gives  50  percent  output  on  the 
SAR  indicator.  Now  there  is  a  difficulty:  the  range  ot  combat  radars  is  a 
closely  guarded  secret  and  to  simulate  noise  properly,  the  range  that  gives  an 
average  noise  output  ot  50  percent  must  be  known.  This  range  can  be  derived 
from  radar  manuals  or  can  be  calculated.  The  simulator  should  also  have 
provisions  to  modify  this  range,  as  the  radar  can  be  operating  at  reduced 
power  or  weather  can  attenuate  the  radar  signal.  Because  the  calculated  range 
can  have  a  slight  error,  adjustments  can  be  easily  made. 

Figure  5.2.1  is  a  block  diagram  ot  a  proposed  noise  generator.  The  look-up 
table  contains  a  normal  distribution  table  that  is  offset  to  give  only 
positive  numbers.  The  hardware  noise  generator  tor  the  counter  is  used  to 
avoid  generating  any  patterns  that  would  occur  it  the  look-up  table  were 
sequenced.  There  are  two  multipliers  that  modify  the  noise  signal.  The  first 
multiplier  modulates  the  noise  signal  as  a  function  ot  range.  The  second 
multiplier  modifies  the  noise  signal  tor  weather,  transmitter  power,  or  to 
change  range  noise  reference.  The  last  function  in  the  noise  generator  is  an 
adder,  where  the  noise  signal  is  added  to  the  simulated  radar  signal. 


5-6 


A  moving  train  attracts  notice  on  the  3.3  nm  patch  map  from  flight  number  59.  Due  to  the  Doppler 
frequency  shift,  the  moving  train  is  displaced  to  the  left  of  the  track.  The  train  Is  moving  along 
the  S-curve  away  from  the  river  at  a  calculated  speed  of  15  mph.  The  drawing  to  the  right  and  the 
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An  example  ot  noise  generator  operation  is  the  best  way  to  understand  this 
concept.  The  parameters  are  a  sector  scan  20  by  40  nautical  miles  with  noise 
equal  to  halt  ot  tull  scale  on  the  display  at  70  nautical  miles. 
Additionally,  heavy  rain  ettects  will  be  calculated  and  will  reduce  the  noise 
range  by  20  nautical  miles  trom  70  nautical  miles  (trom  50  to  70  nm)  as  shown 
in  two  examples: 


Example  one 


Noise  range 


Signal  to  noise  at  20  nm 
Signal  to  noise  at  40  nm 


70  nautical  miles 
150.06 
9.38 


Example  Two 
Noise  range 

Signal  to  noise  at  20  nm 
Signal  to  noise  at  40  nm 


50  nautical  miles 
39.06 
2.44 


The  noise  model  can  be  adjusted  to  any  radar  or  mode  by  changing  the  sottware- 
con trolled  constant  ot  the  second  multiplier. 


5.3  GLITTER 


Synthetic  aperture  radar  and  real  beam  radar  exhibit  a  strong  return 
phenomenon  that  appears  on  one  or  several  scans  and  then  disappears.  This 
phenomenon  is  called  glitter  and  is  caused  when  an  object  on  the  ground 
becomes  a  radar  corner  retlector  tor  a  briet  period.  Figure  5.3.1  is  a 
diagram  ot  an  ideal  three-dimensional  corner  retlector  and  a  two-dimensional 
corner  retlector  and  shows  the  path  ot  the  retlected  waves.  The  three 
surtaces  ot  a  three-dimensional  corner  retlector  need  not  meet.  As  an 
example;  two  buildings  and  the  ground  would  torm  a  three-dimensional  corner 
retlector.  Buildings  separated  by  small  distances  only  reduce  retlectance  by 
a  small  amount.  Three-dimensional  corner  ret  lectors  are  quite  common  when 
mapping  cities  as  most  buildings  are  laid  out  on  grids.  The  return  trom  the 
corner  retlector  in  SAR  is  orders  ot  magnitude  greater  than  in  those  ot  RAR 
because  ot  SAR's  greater  resolution. 
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There  are  two  methods  tor  simulating  glitter.  First,  all  objects  that  cause 
glitter  can  be  added  to  the  data  base.  However,  because  present  data  bases 
are  lacking  in  detail  and  individual  buildings  are  not  coded,  this  method  ot 
glitter  generation  must  await  the  generation  ot  high  resolution  data  bases.  A 
second  method  to  generate  glitter  is  to  randomly  and/or  selectively  put  out  a 
high  level  return,  particularly  in  commercial  and  industrial  areas  where  most 
glitter  is  generated.  In  SAR,  the  glitter  return  will  last  tor  more  than  one 
display.  Since  SAR  is  a  high  resolution  device,  squint  angle  changes  will  be 
minor  in  comparison  to  the  angle  ot  return  tor  the  corner  ret  lector. 

There  is  another  important  aspect  to  glitter  and  high  resolution  SAR.  Study 
Figures  5.3-2,  5.3-3,  5.3-4  and  5.3-5.  Note  that  Figures  5.3-2  and  5.3  3  are 
photographs  ot  Busch  Stadium  in  St.  Louis,  Missouri.  The  product 
specit icat ion  tor  the  Digital  Landmass  System  (DLMS )  data  base  has  a  code 
(321)  tor  an  enclosed  stadium.  At  the  present  time,  most  radar  simulators 
generate  a  point  tor  DLMS  code  321.  For  SAR,  code  321  should  be  decoded  as  a 
ring  with  a  vertical  outside  wall,  a  sloping  interior,  and  glitter  returns  at 
the  locations  ot  light  towers.  Figure  5.3-4  is  a  SAR  display  ot  the  Daggett 
Marine  Corps  depot.  The  DLMS  code  tor  warehouses  is  861.  For  Daggett  Marine 
Corps  warehouses,  SAR  simulation  requires  insertion  ot  glitter  warehouse  roots 
where  vents  and  air  conditioning  units  are  located.  The  DLMS  code  that 
generates  glitter  should  be  identitied  tor  proper  SAR  simulation.  One  rule  ot 
thumb  tor  surtace  codes  in  the  DLMS  codes  is  that  surtaces  which  are  highly 
retlective  generate  weak  returns  or  generate  glitter.  The  codes  can  also  be 
omniret lective.  DLMS  code  155  is  a  heavy  tabrication  industrial  building  with 
a  sawtooth  root.  In  low  grazing  angles,  this  root  will  retlect  trom  steep 
sides  ot  the  sawtooth. 

5.4  SCINTILLATION  AND  MULTILOOK 

SAR  has  another  ettect  that  ditters  trom  conventional  pulse  radar  known  as 
scintillation  or  speckle.  SAR,  as  a  coherent  system,  is  very  sensitive  to  any 
change  in  the  transmission  path  length  and  target  movement  and  multipath 
returns  also  downgrade  SAR  pertormance.  The  quality  ot  the  atmosphere  is  not 
homogeneous,  but  variable  and  subject  to  constant  changes.  The  twinkling  ot 
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stars  is  an  example  ot  scintillation.  Some  SAR  systems  reduce  scintillation 
by  a  tradeott  with  resolution.  This  tradeott  is  accomplished  by  averaging  a 
number  ot  returns  ot  the  same  area.  For  combat  SAR's  the  number  used  tor 
averaging  is  generally  two  or  tour. 

Nonmobile  target  motion  is  caused  by  wind.  Strongly  built  structures  do  not 
give  any  notable  scintillation  but  vegetation  driven  by  wind  can  cause 
scintillation,  as  well  as  smearing  the  return.  As  would  be  expected,  grass 
generates  more  scintillation  than  trees;  water  varies  trom  nonret lective  to  a 
highly  reflective  return. 

The  simulation  ot  wind,  rough  water  scintillation,  and  multilook  are 
implemented  in  the  same  manner.  Simulation  begins  with  an  unmodified 
display.  The  overall  picture  can  then  be  modified  by  using  a  Gaussian 
multiplier  look-up  table  to  generate  atmospheric  scintillation;  there  should 
be  three  look-up  tables  (i.e.,  tor  one,  two,  or  tour  looks).  It  rough  water 
is  identifiable  in  the  data  base,  a  random  look-up  table  can  be  used  to  vary 
the  returns  and  the  returns  should  also  be  smeared.  A  coastline  return  can  be 
modified  by  rough  water  in  an  area  such  as  the  Columbia  River  Bar.  Vegetation 
driven  by  wind  is  simulated  by  modifying  returns  with  a  look-up  table  and 
smearing.  For  high  or  gusty  winds  all  the  look-up  tables  could  be  used.  Low 
velocity  winds  can  be  simulated  by  not  modifying  data  on  some  scans. 

Figure  5.4-1  is  a  photograph  ot  two  SAR  displays.  The  first  display  is  made 
up  ot  six  subscans  and  has  a  mult  Hook  ot  one.  Note  that  fading  is  due  to  the 
antenna  footprint  and  the  amount  of  glitter.  The  second  display  is  ot  the 
same  area  except  that  the  display  has  a  multilook  ot  two.  Fading  and  glitter 
have  been  reduced.  simulating  multilook  tor  the  antenna  footprint  is  not 
complicated.  For  Multilook  1  use  the  antenna  footprint  multiplication  table 
calculated  from  the  actual  radar  to  modify  picture  returns;  for  Multilook  2 
use  a  modified  antenna  footprint  table  to  modify  picture  returns;  and  tor 
Multilook  4  do  not  modify  the  picture. 
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5.5  FOOTPRINT 

The  area  size  that  can  be  processed  by  SAR  is  limited  to  the  area  size 
illuminated  by  the  antenna.  For  SAR  processing,  this  is  generally  considered 
to  be  the  area  within  the  two-way  path,  3  dB  points. 

There  are  generally  two  types  ot  antenna  used  in  aircratt  SAR.  For  patch 
mode,  a  pencil  beam  is  considered  best  as  it  concentrates  the  most  RF  energy 
in  the  area  ot  Interest.  The  sector  scan  mode  requires  a  long  and  narrow  area 
to  be  illuminated.  A  tan-shaped  beam  is  used  tor  this  purpose.  Figure  5.5-1 
is  a  diagram  ot  an  aircratt  antenna  tootprint. 

Simulation  ot  the  antenna  tootprint  is  relatively  simple.  The  antenna 
tootprint  algorithm  consists  ot  multiplying  the  reflectance  values  by  a 
coetticient  that  is  a  tunction  ot  the  azimuth  angle.  A  simple  approach  is  to 
assume  that  the  coetticlents  are  constant  with  slant  range.  This  is  not 
entirely  true,  but  tor  simulation  purposes  it  is  a  valid  assumption.  Figure 
5.5-2  is  a  diagram  ot  a  typical  antenna  pattern.  Using  Figure  5.5-2,  a 
velocity  error  example  can  be  demonstrated.  Figure  5.5-3  shows  the  amplitude 
response  with  velocity  errors.  It  signiticantly  great  velocity  error  occurs, 
the  display  could  become  blank. 

The  velocity  error  algorithm  consists  ot  a  multiplier  and  a  look-up  table. 
The  pointer  to  the  look-up  table  is  advanced  trom  scanline  to  scanline  and  is 
ottset  by  velocity  errors.  Coetticlents  that  are  in  the  look-up  table  will  be 
determined  by  the  aircratt' s  antenna  and  the  mode  ot  the  SAR. 

5.6  DOPPLER  ELEVATION  MAPPING 

In  previous  sections,  it  has  been  stated  that  SAR  mapping  is  pertormed  on 
isodoppler  lines  and  that  the  theory  ot  SAR  developed  by  considering  that  the 
earth's  plane  is  tlat.  We  now  look  at  the  problem  that  occurs  when  points  ot 
interest  are  not  in  the  reference  plane,  such  as  tall  buildings,  mountains  and 
depressions.  The  information  that  is  needed  when  considering  SAR  elevation 
mapping  is  a  determination  ot  how  the  aircratt  SAR  determines  elevation  ot  the 
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SAR  reterence  plane.  The  relative  velocity  ot  the  center  point  ot  the  scan  on 
the  ground  and  the  aircratt  is  used  as  a  reterence  in  aircratt  SAR 

processing.  The  algorithm  that  determines  this  velocity  should  be  a  guide 

when  determining  what  algorithm  will  determine  the  elevation  reterence  plane. 
One  type  ot  algorithm  is  to  average  a  number  ot  points.  For  SAR  elevation 
mapping,  we  will  make  an  assumption  that  the  elevation  reterence  plane  is 
known. 

A  comparison  has  been  made  ot  SAR  to  a  large  lens  or  antenna.  For  all 
calculations,  the  center  ot  this  lens/antenna  is  considered  to  be  haltway 
along  the  data  collection  path;  the  nadir  point  is  the  origin  ot  the 

isodoppler  map.  Now  each  point  in  the  isodoppler  map  has  its  own  signature  in 
azimuth  depression  angle  and  slant  range.  The  isodoppler  lines  are 

represented  by  Equation  5.6.1. 


Lid  =co6(\(i)  co6(0w) 


Where 

L  =  an  isodoppler  line 


(5.6.1) 


The  depression  angle  is  represented  by  Equation  5.6.2: 
[ Arp  -  Erp 


V’  =  taiT 


(5.6.2) 


Where 

A 

RP 

E 

RP 


=  altitude  above  the  reterence  plane 
=  elevation  above  the  reterence  plan 
=  ground  range 


Let  the  depression  angle  ot  the  ground  range  be  the  reterence  when  the  value 
ot  elevation  above  the  reterence  plane  is  zero.  Since  SAR  processing  is 
pertormed  along  isodoppler  lines,  the  equation  tor  a  point  that  is  not  in  the 
reterence  plane  is: 
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(5.6.3) 


Where 


6 

AM 


«  mapping  azimuth  trom  the  velocity  vector 


Solving  tor  the  azimuth  angle  0  gives  the  azimuth  location  of  a  point 

AM 

that  is  not  at  the  elevation  ot  the  reference  plane. 


7  AM 
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os(0  w) 


cos(^) 
cos(  V>fi ) 


(5.6.5) 


To  indicate  how  Equation  5.6.5  varies  with  azimuth  and  depression  angle. 
Equation  5.6.5  can  be  differentiated. 

dOjMp  Un(V’fl) 

-  - - —  \  J  .D  .0  J 

d  V'/j  d  0  ^  tan(  0  w ) 


This  type  ot  equation  has  been  shown  before;  it  has  the  same  form  as  the  rate 
ot  change  ot  the  isodoppler  lines,  as  one  would  expect.  An  object  that  is 
above  the  reference  plane  will  plot  in  azimuth  toward  the  velocity  vector  and 
an  object  that  is  below  the  reference  plane  will  plot  away  from  the  velocity 
vector.  Figure  5.6-1  is  a  plot  ot  the  amount  ot  shift  in  azimuth  tor  a 
1000- toot  object,  at  10  nautical  miles. 

This  layover  is  peculiar  to  SAR  and  has  several  implications  in  SAR 
simulation.  One  characteristic  is  that  the  shadow  ot  an  object  with  layover 
stays  where  it  is.  Another  characteristic  is  that  an  object  with  layover  maps 
into  other  objects  or  shadows.  Figure  5.6-2  is  a  diagram  ot  a  1000-foot 
object  and  an  indicator  ot  the  layover  phenomenon.  The  difference  in  slant 


A  A 


NAVTRAEQUIPCENT  83-C-0068 


range  between  the  top  and  bottom  ot  the  building  is  251  teet.  The  top  ot  the 
object  maps  toward  the  aircratt  and  also  toward  the  velocity  vector.  The 
layover  mapping  ot  terrain  (mountains)  requires  only  one  pixel  to  be  mapped. 
A  tall  object  must  be  represented  by  a  number  ot  pixels  and  each  pixel  must  be 
remapped.  Layover  can  be  simulated  but  this  requires  a  considerable  amount  ot 
calculation. 


The  specification  ot  the  layover  function  in  a  SAR  simulator  will  depend  upon 
three  factors:  1)  the  degree  ot  realism  that  is  required,  2)  structural 
heights  ot  the  area  to  be  simulated,  and  3)  SAR  resolution.  Two  examples  are 
included  as  a  basis  tor  the  decision-making  process: 

Example  One 


Range 

Pixel  size 
Resolution 
Height  ot  structure 
Azimuth 

Depression  angle 
Layover 

Example  Two 
Range 

Pixel  size 
Resolution 
Height  ot  structure 
Azimuth 

Depression  angle 
Layover 


10  nautical  miles 
10  teet 
20  teet 
1200  teet 
15  degrees 
15  degrees 
1.06  degrees 
112  pixels 

10  nautical  miles 
50  teet 
100  teet 
150  teet 
15  degrees 
15  degrees 
0.13  degrees 
3  pixels 
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In  the  tirst  example;  without  layover  the  simulator  display  would  be  a  poor 
match  to  the  real  world  display.  In  the  second  example,  the  simulator  could 
produce  a  good  match  without  layover  simulation.  A  word  ot  caution  should  be 
inserted  here:  all  modern  radars  are  designed  to  be  upgraded  and  even  it 
layover  is  not  specitied,  a  capability  tor  inserting  overlay  at  a  later  date 
should  be  designated. 


5.7  MALFUNCTION  EFFECTS 


For  proper  simulation  velocity  errors  must  be  included.  In  the  principles  ot 
SAR  it  war.  emphasized  that  velocity  was  a  linear  phase  shitt  across  the 
aperture.  This  linear  phase  shitt  causes  the  aperture  to  be  steered.  Figure 
5.7-1  and  5.7-2  are  examples  ot  steering  the  aperture  out  ot  the  area  being 
illuminated  by  the  radar  antenna.  Figure  5.7-3  is  a  diagram  indicating 
velocity  error.  A  solution  tor  velocity  error  simulation  is  given  in  Equation 
5.7.1. 


Fds  —  ■ 


2  VacF 


(5.7.1) 


The  solution  tor  azimuth  shitt  is  to  equate  two  Doppler  trequencles  that  have 


dltterent  aircratt  velocities: 


cos(0kvt)cos(V’)  K1c  =co8(8Kv  )cos(  V'H  vac  + 


(5.7.2) 


Where 


GvvE  =  az*muth  trom  the  velocity  vector  when  a  velocity  error  exists 

v  =  velocity  error 

ERR 


Solving  tor  the  azimuth  angle  tor  a  velocity  error  gives  Equation  5.7.3. 


0  WE  =*COS  *  |cOs(  0  w) 


(  K/tc  +  Verp  ) 


(5.7.3) 
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A  velocity  used  tor  SAR  calculations  that  is  higher  than  the  actual  velocity 
will  shilt  the  aperture  toward  the  velocity  vector.  A  velocity  error  on  the 
low  side  will  shitt  the  aperture  away  from  the  velocity  vector.  An  example  ot 
velocity  error  is  provided: 

Example:  Velocity  Error 


Azimuth  angle 
Range 
Velocity 
Velocity  Error 


=  20  degrees 

=  10  nautical  miles 

=  300  knots 
=  10  knots 


This  10  knot  velocity  error  shifts  the  aperture  6.17  degrees  toward  the 
velocity  vector.  At  a  range  ot  10  nautical  miles  this  is  a  shitt  ot  6570 
teet,  emphasizing  that  a  small  velocity  error  makes  a  large  SAR  processing 
error. 


To  tind  out  how  azimuth  error  varies  as  a  function  ot  the  azimuth  angle  and 
velocity,  Equation  5.7.4  is  ditterentiated  with  respect  to  the  velocity  error: 

^  _ ! _  (5.7.4) 

Werr  Vxc'tan(©v'v) 


Equation  5.7.4  shows  that  velocity  errors  give  a  large  azimuth  error  near  the 
velocity  vector. 

The  simulation  ot  velocity  errors  is  straightforward.  First,  consider  the 
center  ot  the  radar  picture  as  a  reference.  The  azimuth  shitt  in  this 
reference  pixel  due  to  velocity  error  is  now  calculated.  This  azimuth  shitt 
is  now  used  to  calculate  a  set  ot  new  data  retrieval  coordinates.  This  new 
data  is  then  used  to  generate  a  radar  display.  Figure  5.7-1  is  an 
illustration  ot  data  shitt  due  to  a  velocity  error. 

Velocity  error  causes  a  lading  or  shading  ot  the  radar  picture.  During  normal 
operation  ot  SAR,  the  area  ot  interest  occupied  is  within  the  3  dB  points  ot 
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the  antenna  tootprint.  A  velocity  error  shltts  the  designed  display  area  out 
ot  the  display  and  brings  in  an  area  that  is  outside  the  3  dB  points  ot  the 
antenna  tootprint.  In  the  simulator  the  solution  tor  tading  due  to  velocity 
error  Is  to  shift  the  antenna  tootprint  multiplying  table  to  match  the 
velocity  error. 

Acceleration  is  detined  as  the  derivative  ot  velocity  relative  to  time. 
Acceleration  errors  can  be  caused  by  three  related  factors.  The  tirst  t actor 
could  be  a  taulty  inertial  system;  the  second  tactor  could  be  a  maneuver  ot 
the  aircratt  that  exceeds  the  inertial  system's  correction  capability.  The 
third  tactor  is  tlexing  ot  the  aircratt  that  causes  motion  between  the 
inertial  system  and  the  radar  antenna. 


Acceleration  errors  cause  loss  ot  resolution  in  the  azimuth  direction. 
Acceleration  simulation  is  straightforward.  This  simulation  consists  ot 
spreading  the  retlectance  data  ot  one  pixel  over  a  number  of  pixel  locations. 
A  small  acceleration  error  would  involve  only  a  tew  pixel  locations  and  a 
large  acceleration  error  would  involve  many  pixels.  Most  acceleration  errors 
in  an  actual  aircratt  are  very  complex  and  would  be  Impossible  to  simulate. 
To  understand  acceleration  errors,  a  linear  acceleration  can  be  used. 


I 

The  starting  point  will  be  Equation  5.7.1  (a  velocity  error  equation).  This  i 

equation  is  rewritten  to  include  acceleration  errors; 


Fd 


+  /  /Mcc) 


(5.7.5) 


For  example,  the  acceleration  function  is  set  to  a  constant.  This  is  a  linear 
acceleration  error.  Equation  5.7.6  is  the  equation  tor  a  linear  acceleration 
error; 

Fd  ’“2  — |  VAC  +  Accih*  ^l)! 


5-33 


•-W- 


v.w*rv 


'  .  *  .  * 
*  «  *  ,  •  •  ' 


V*V*  .*•  »**  v  * 

.*  *  •  o  o  •  • .  •  * » •  o  •:  v 


(5.7.6) 


i 


NAVTRAEQU 1PCENT  83-C-0068 


At  the  time  t^  the  Doppler  frequency  is  2*F»Vftc/C  and  at  t^  the 


Doppler  frequency  is  2*F/CLV  +A  (t  -t  )J. 

AC  CC  2  1 


By  replacing  the  velocity  error  in  Equation  5.7.1  with  the  acceleration  error, 
the  equation  for  the  azimuth  angle  of  the  end  of  the  data  collection  can  be 
obtained. 


0  WF  =COS  1 


cos(0kv)  |  K«C  +  Acc(h  ~  <l)| 


1  AC 


(5.7.7) 


The  azimuth  angle  at  the  start  of  the  data  collection  is  8  .  The  data 

vv 


tor  one  pixel  is  spread  over  the  azimuth  from  0  tc  0  . 

vv  WE 


An  example  of  the  effects  of  acceleration  errors  is  given: 


Example:  Acceleration  Error 


Azimuth  angle 

Range 

Velocity 

Acceleration 

Start  of  data  collection 

End  of  data  collection 

Azimuth  angle  at  the 

end  of  data  collection 

Pixel  spread 


20  degrees 
10  nautical  miles 
300  knots 

1.2  knots  per  second 
0  seconds 
10  seconds 


12.23  degrees 
7.76  degrees 


Although  this  example  is  quite  drastic,  it  does  illustrate  what  acceleration 
errors  can  do  to  a  SAR  image.  In  this  case,  a  real  beam  radar  would  give  a 
better  picture. 
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6.0  DATA  BASE  REQUIREMENTS  FOR  SAR  SIMULATION 

6.1  ELEVATION  DATA  BASE 

Requirements  tor  SAR  elevation  data  are  more  exacting  than  those  tor  real 
aperture  radar  because  SAR  resolution  is  trom  16  to  more  than  100  times 
greater.  For  this  reason,  short  cuts  that  work  tor  RAR  tall  short  in  SAR. 
The  Digital  Landmass  System  (DLMS)  data  base  tor  elevation  specities  elevation 
data  on  3  second  by  3- second  posts  trom  0  to  50  degree  latitude.  At  the 
equator  this  would  be  303.6  teet  by  303.6  teet.  A  high  resolution  SAR 
operating  at  10- toot,  resolution  with  a  display  ot  230,400  pixels  uses  only  256 
posts  to  generate  all  ot  the  elevation  data.  The  above  tigures  were 
calculated  tor  a  location  ot  zero  latitude.  Fortunately  there  are  a  number  ot 
things  that  can  be  done  to  make  the  DLMS  usetul  tor  SAR. 

The  most  simple  way  to  retine  the  DLMS  elevation  data  base  is  to  use  a  two 
dimensional  linear  interpolation  algorithm  to  obtain  elevation  data  tor  points 
between  posts.  Figure  6.1-1  depicts  the  two-dimensional  linear  interpolation 
algorithm.  However,  this  algorithm  has  two  shortcomings:  1)  the  elevation 
slope  is  essentially  made  up  ot  linear  lines,  and  2)  the  elevation  slope  has 
discontinuities.  These  discrepancies  can  be  tolerated  in  some  RAR  systems  but 
are  not  acceptable  in  high  resolution  SAR  systems.  To  remove  the  straight 
lines  and  discontinuities,  the  derivative  ot  the  slope  is  required  to  be  a 
continuous  tunction.  The  use  ot  16  posts  and  the  two-dimensional  weighted 
parabolic  interpolation  algorithm  does  this.  There  are  other  algorithms  that 
would  do  the  same  tunction,  but  the  two-dimensional  weighted  parabolic 
interpolation  algorithm  is  one  ot  the  most  simple.  Figure  6.1-2  is  the 
diagram  and  equations  ot  the  two-dimensional  weighted  parabolic  interpolation 
algorithm.  This  algorithm  and  any  other  algorithms  will  give  an  overshoot  tor 
a  step  tunction.  Figure  6.1-3  is  an  example  ot  the  overshoot.  The  maximum 
overshoot  tor  a  two-dimensional  weighted  parabolic  interpolation  algorithm  is 
0.0741  ot  the  step.  The  DLMS  data  base  contains  very  tew  step  tunctions, 
hence  the  step  tunction  response  should  not  be  a  problem. 


EA  =  0.5  [ (-E1+3E2-3E3+E4)y33 
+0.5  t(2E1-5E2+4E3-E4)y23 
+0.5  C (E3-E1 )  y  3+E2 
EB  =  f  [E5,E6,E7,E83y 
EC  =  f  [E9.E10.E1 1 ,E123  y 
ED  =  f  [E13.E14, E15.E163 y 
E  =  f  [EA,EB,EC, ED3X 


Figure  6.1-2  PARABOLIC  INTERPOLATION  ALGORITHM 
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High  resolution  SAR  has  tiner  shadow  details  than  can  be  obtained  by  just 
interpolating  the  DLMS  data  base.  Figures  6.1-4  and  6.1-5  show  how  much 
shadow  detail  there  is  in  low  grazing  angle  high  resolution  SAR.  There  is  a 
simple  solution  to  make  the  elevation  data  base  more  detailed  with  little  loss 
ot  accuracy  at  the  data  posts.  Only  at  the  data  posts  is  the  elevation  known 

and  this  data  has  a  tolerance  assigned  to  it.  By  putting  in  a  roughness 

algorithm  that  modities  the  elevation  trom  the  two-dimensional  weighted 
parabolic  interpolation  algorithm,  the  tine  variations  in  the  elevation  data 
can  be  obtained  with  little  loss  in  accuracy.  Figure  6.1-6  is  a  block  diagram 
ot  how  a  roughness  algorithm  can  be  implemented.  A  random  delta  X  generator 
and  a  random  delta  Y  generator  are  the  starting  points  tor  this  algorithm; 

magnitude  and  time  ot  change  are  both  random.  Each  delta  is  put  into  an 

integrator  and  the  outputs  are  summed.  The  sum  is  subject  to  an  out- ot- limits 
check.  It  the  sum  is  out  ot  limits,  the  sum  is  truncated  to  tfte  limit.  The 
output  ot  the  limiter  Is  now  summed  with  the  elevation.  This  algorithm  can  be 
made  coherent  with  the  terrain  location. 

At  this  time  the  DLMS  elevation  data  base  does  not  provide  tor  a  roughness 
coet tic  lent  In  terrain.  A  clue  to  roughness  can  be  obtained  trom  the 
elevation  data  base.  The  amount  ot  elevation  moditicatlon  is  a  tunction  ot 
the  terrain.  For  example,  take  three  locations  and  assign  a  roughness  index. 
The  tirst  location  considered  Is  wheat  tields  ot  Kansas.  In  this  case  the 

roughness  index  will  be  zero,  hence  the  elevation  will  not  be  moditied.  The 

second  location  considered  Is  toothills  ot  the  Coast  Range  Mountains  ot 
California.  This  location  should  be  assigned  a  medium  roughness  index.  The 
last  example  is  Bad  Lands  ot  South  Dakota.  The  roughness  index  tor  this 
location  should  be  a  high  index.  Although  the  elevation  that  is  generated  by 
the  roughness  algorithm  is  not  an  exact  reproduction  ot  the  actual  terrain, 
the  added  realism  gained  by  using  the  roughness  algorithm  more  than  makes  up 

tor  any  loss  ot  accuracy  in  elevation  at  the  posts.  Figure  6.1-7  is  an 

example  ot  the  roughness  algorithm. 
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Figure  6.1-7  EXAMPLE  OF  ROUGHNESS  ALGORITHM 
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6.2  PLANIMETRIC  DATA  BASE 

The  best  algorithms  tor  radar  simulation  will  not  give  a  good  simulated  radar 
display  unless  these  algorithms  have  a  realistic  planimetric  data  base  tor 
computation.  For  this  reason,  the  planimetric  data  base  is  the  heart  ot  radar 
simulation. 

In  1972  the  Detense  Mapping  Agency  (DMA)  was  directed  to  produce  a  digital 
data  base  tor  radar  simulation.  In  this  same  time  trame.  Link  built  a  digital 
radar  simulator  tor  the  U.S.  Air  Force.  This  simulator  (the  1183)  used  the 
DMA  data  base  directly.  In  1974  DMA  published  the  production  specit ication 
(DLMS)  tor  digital  data  base  generation  ot  radar  data.  The  1183  simulator  and 
the  data  base  generated  by  DMA  were  tested  and  evaluated  by  the  Air  Force, 
DMA,  and  Link.  The  DLMS  specit  ication  was  revised  in  July  ot  1977  and  again 
in  June  ot  1983.  In  1978  the  Air  Force  requested  DMA  to  generate  a  data  base 
to  support  visual,  intrared  (IR),  and  radar.  This  request  resulted  in 
publication  ot  a  product  specit ication  tor  a  prototype  data  base  to  support 
high  resolution  sensor  simulation  in  December  ot  1979.  The  B-lB  SAR 
requirements  have  been  responsible  tor  generation  ot  a  product  specit ication 
to  support  high  resolution  Level  X  data  bases  tor  synthetic  aperture  radar, 
ettective  in  June  ot  1983.  DMA  has  developed  two  areas  ot  Level  X  data  bases 
tor  evaluation  and  is  currently  working  on  a  third.  At  this  time,  DMA  has 
generated  a  moditied  Level  I  test  data  base  tor  evaluation  ot  SAR  simulation. 
The  specit ication  title  is:  DMA  Product  Specit ication  Supporting  Prototype 
Level  X  tor  STRC,  January  1984,  Stock  No.  Spec  X  HRD/B1B.  The  DMA  high 
resolution  ettort  is  being  conducted  tor  the  Air  Force  in  support  ot  the  B-lB 
radar  simulator. 

A  comparison  ot  Level  I,  Level  II,  and  Level  X  data  bases  is  given  to  Indicate 
the  degree  ot  data  base  support  tor  SAR  simulation.  Minimum  limits  tor  coding 
are  listed: 
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CATEGORY :  ROADS 

.1 

i 

Lenqth 

Width 

Heiqht 

rV*V 

Level  I 

300m 

Any 

Any 

•*. 

.  '  »  »  1 
-  •  •  • 

Level  II 

300m 

Any 

Any 

“  a  *■». 

Level  X 

10m 

Any 

Any 

r? 

K  ... 

CATEGORY:  BRIDGES 

Lenqth 

Width 

Heiqht 

Level  I 

30m 

Any 

Any 

-  -1 

Level  II 

30m 

Any 

Any 

Level  X 

10m 

Any 

Any 

•Y-Y 

CATEGORY :  ISOLATED  STRUCTURES 

Lenqth 

Width 

Heiqht 

Level  I 

30m 

Any 

3m 

d 

Level  II 

10m 

Any 

3m 

Level  X 

10m 

Any 

There  is  not  a  great  ditterence  between  Level  I  and  Level  II.  Level  X, 
however,  is  ot  much  greater  resolution  than  Level  I  or  Level  II.  At  this 
time,  the  Detense  Mapping  Agency  and  the  Air  Force  are  trying  to  resolve  the 
requirements  tor  SAR.  The  amount  ot  Level  X  data  will  be  limited  because  ot 
the  high  cost  ot  generation.  Where  there  is  a  lack  ot  Level  X  data,  Level  I 
or  Level  II  data  will  have  to  suttice.  Only  with  Level  X  data  will  a  SAR 
simulated  image  be  ground  truth.  DMA  will  be  the  source  tor  all  production 
data,  but  a  short  turnaround  tor  a  test  or  quick  update  and  small  areas  ot 
data  bases  can  be  generated  on-site.  The  problem  ot  enhancing  Level  I  data  to 
give  a  realistic  SAR  image  is  a  large  one.  These  realistic  images  will  not  be 
ground  truth. 


Vl 
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A  short  explanation  ot  some  terms  used  in  radar  simulation  is  included  prior 
to  discussion  ot  ways  to  enhance  the  data  base.  The  tirst  two  radar  simulation 
terms  to  be  discussed  are  retlectance  and  signature.  In  the  past,  only  a 
retlectance  was  calculated  and  the  signature  was  Ignored.  An  example  is  given: 
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EXAMPLE:  LEVEL  I  CODE 


Feature  analysis  N/A 
Feature  type  (aerial)  2 
Predominant  height  (24  feet)  4 
Feature  identification  code  (houses)  420 
Surface  material  category  code  (composition)  4 
Number  of  structures  per  square  kilometer  (680-825)  5 
Percent  tree  cover  (20%)  2 
Percent  root  cover  (20%)  2 
An  outline  ot  the  area 


Present  simulators  generate  a  reflectance  code  based  on  the  surface  material 
category.  DMA  has  provided  all  the  information  necessary  to  generate  a 
realistic  signature  ot  a  residential  neighborhood  with  trees  and  streets.  The 
signature  would  be  realistic  but  it  would  not  be  ground  truth.  This  realistic 
signature  could  be  generated  by  at  least  four  methods.  The  first  method 
expands  the  DMA  data  base  oft  line.  The  second  method  expands  the  data  base 
on  line.  The  third  method  uses  an  interim  butter  and  interlocking  overlays  to 
generate  the  signatures  with  hardware.  The  residential  neighborhood  signature 
has  tour  retlectance  levels.  These  levels  are  ground,  trees,  rooftops,  and 
roads.  A  RAR  simulator  makes  use  Of  one  level  of  retlectance  coding  tor 
medium  and  long  ranges.  But  at  short  ranges,  RAR  simulators  lack  realism  due 
to  the  lack  ot  signatures.  For  SAR  simulation  a  good  signature  is  required. 
Every  feature  identification  number  has  its  own  signature  that  consists  ot 
more  than  one  retlectance  level. 

The  predominant  height  and  individual  height  have  been  used  interchangeably  in 
the  past.  To  get  good  SAR  simulation,  the  predominant  height  must  be  broken 
up  in  a  random  way  as  a  function  of  the  feature  identification  number.  A 
residential  tract  would  have  very  little  modification  ot  predominant  height, 
but  an  area  ot  commercial  buildings  would  have  large  variations  in  predominant 
height.  Much  work  needs  to  be  done  in  the  expansion  ot  DMA  feature  codes. 
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6.2.1  China  Lake  Air  Base  Data  Base  Generation 

The  generation  ot  a  high  resolution  data  base  tor  China  Lake  Air  Base  was 
straight  forward .  The  first  task  was  to  generate  a  computer  program  that  would 
allow  input  ot  data  through  a  console  into  a  tile  that  matched  the  Link  32-bit 
word,  Level  I  format.  The  next  task  was  to  obtain  all  information  that  was 
available  on  the  China  Lake  Air  Base  region.  The  following  information  was 
obtained  from  the  U.S.  Geological  Survey  in  Menlo  Park,  California. 

1)  Ridgecrest  North,  California:  map  scale  1/24000 

2)  Ridgecrest  North,  California:  orthophotograph 

3)  Ridgecrest  aerial  photograph 

The  area  around  the  airport  was  enlarged  photographically  by  a  factor  ot 
tour.  The  enlargement  scale  factor  is  1/6000. 

This  means  that  one  inch  on  the  enlargement  is  equal  to  500  feet.  The  program 
contained  a  scale  factor  adjustment  that  was  accurate  to  one  part  in  1/100000, 
but  it  was  not  required  as  the  enlargement  error  was  less  than  one  part  in 
10000.  A  gi idded  vellum  was  placed  on  the  enlarged  photograph.  The  lower 
lett  corner  ot  the  vellum  was  referenced  to  the  DMA  manuscript  that  covers  the 
China  Lake  Air  Base.  The  major  task  was  to  identity  the  targets  on  the  map 
and  code  them  according  to  the  DMA  Level  I  product  specification.  The  output 
ot  the  tile  that  was  manually  coded  was  plotted  out  and  the  errors  were 
Identified  and  corrected.  The  DMA  manuscript  was  searched  and  the  features 
that  conflicted  with  the  manually  generated  data  base  were  deleted.  The  last 
step  was  to  merge  the  manually  generated  data  base  with  the  DMA  data  base. 
Figure  6. 2. 1-1  is  an  F- 15  SAR  picture  ot  the  China  Lake  Air  Base. 


Photograph  of  China  Lake 
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6.3  COMMENTS  ON  SAR  PHOTOGRAPHS 

i 

i 

6.3.1  Comments  on  Hard/Sot t  Target  Switch 

Reterence:  Figure  3.8.4 

I 

These  two  photographs  are  a  good  example  at  what  happens  when  SAR  saturates. 
The  azimuth  is  smeared  across  one- third  ot  the  width  of  the  picture.  The 

smearing  stops  abruptly  which  would  indicate  that  the  total  picture  consists 
8  ot  three  processed  SAR  pictures  that  were  merged  together.  The  azimuth 

resolution  is  determined  by  SAR  processing  and  the  range  resolution  is 
determined  by  pulse  compression.  These  two  processes  do  not  show  the 

saturation  characteristics.  Pulse  compression  saturation  shows  a  slight 
amount  ot  smear  around  the  target  and  detinite  side  lobes.  A  good  guess  would 
put  the  number  ot  pixels  in  the  range  direction  at  320.  The  saturation  target 
is  in  the  middle  ot  the  picture.  In  the  Hughes  Reports,  the  pulse  compression 
ratio  has  been  ret  erred  to  as  169.  Halt  ot  320  is  160;  hence  it  would  be 
I  expected  that  the  range  saturation  would  extend  across  the  picture.  The 

saturation  does  extend  the  full  range  length  ot  the  picture. 

6.3.2  Comments  on  Low  Squint  Angle 

i 

Reterence:  Figure  4.1.1 

The  three  SAR  pictures  in  this  tigure  illustrate  that  a  modern  SAR  can  map 
'  within  eight  degrees  ot  the  velocity  vector  to  50  degrees  trom  the  velocity 

vector. 

6.3.3  Comments  On  Moving  Ground  Target  (Mudd’s  Landing) 

Reterence:  Figures  5.1.2  and  5.1.3 

This  radar  photograph  is  ot  two  trains.  One  train  is  stationary  and  the  other 
is  moving.  The  stationary  train  has  a  very  strong  return  with  a  radar  shadow 
directly  behind  it.  The  second  train  is  moving  and  leaves  a  shadow  at  the 
same  place  it  would  be  it  it  were  stationary,  except  that  the  shadow  is  longer 
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and  tapered  due  to  movement  along  the  track.  Note  the  return  extends  out  ot 
the  picture  whiLe  the  shadow  is  within  the  picture.  Uneven  movement  breaks  up 
the  return  and  the  return  is  stretched  out  and  dimmed  where  the  train  is  going 
around  a  curve.  All  ot  the  above  mentioned  teatures  can  be  simulated,  but  the 
simulation  can  be  complex. 

6.3.4  Comments  On  0.67  Nml  Patch  Map  Ot  Busch  Stadium 

Reterence:  Figures  5.3.2  and  5.3.3 

These  two  tigures  are  a  challenge  to  the  data  base  designer.  The  stadium 
light  towers  seem  to  show  and  the  press  box  is  a  good  radar  ret  lector.  There 
is  a  great  deal  ot  ditterence  in  retlectivity  among  the  buildings.  The  arch, 
it  it  did  retlect  any  energy,  is  not  recognizable.  Note  the  Stoutter  Hotel 
consists  ot  a  number  ot  bright  ret  lectors  and  they  map  in  several  range 
buckets.  There  seem  to  be  several  pixel  displacements  in  the  Stoutter  Hotel 
due  to  layover. 

6.3.5  Comments  On  Daggett  Marine  Corps  Supply  Center  Solar  One  Power  Plane 

Reterence:  Figures  5.3.4  and  5.3.5 

These  two  SAR  photographs  illustrate  how  large  buildings  become  complex 
targets  at  high  resolution.  A  close  look  shows  the  roots  to  be  gabled  roots 
with  a  monitor.  These  roots  can  be  corner  ret  lectors.  There  also  seems  to  be 
a  large  number  ot  root  vents  that  are  strong  ret  lectors.  In  the  17-toot 
resolution  SAR  picture,  the  long  building  that  taces  the  aircratt  is  such  a 
strong  retlector  that  it  is  saturating  the  SAR.  In  tact,  quite  a  number  ot 
targets  are  saturating  the  SAR.  There  seems  to  be  a  chain  link  tence  around 
the  base  that  is  a  strong  retlector. 
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6.3.6  Comments  On  Radar  Map  Quality 

Reterence:  Figure  5.4.1 

These  two  photographs  show  what  effect  the  multilook  algorithm  has  on  picture 
quality.  The  photographs  are  composed  of  six  subpictures.  A  multilook  ot  one 
gives  a  scalloping  ettect  (Venetian  blind  effect)  due  to  the  antenna 
footprint.  When  two  pictures  are  combined  and  the  footprints  are  staggered, 
the  scalloping  is  eliminated.  This  is  a  multilook  ot  two;  multilooks  greater 
than  two  have  been  generated.  The  simulation  world  differs  from  the  real 
world  in  that  you  start  with  a  good  quality  picture  and  put  the  scalloping  in. 

6.3.7  Comments  On  PVU  Correction  Ot  Antenna  Shading 

Reterence:  Figures  5.7.1  and  5.7.2 

This  radar  picture  is  an  excellent  example  of  velocity  error  in  SAR.  The 
simulation  of  the  phenomenon  ot  SAR  should  be  a  requirement  as  it  will  train 
the  pilot  to  update  the  altitude  that  is  used  in  SAR  processing,  as  well  as 
checking  on  the  inertial  platform.  The  velocity  error  has  been  covered  in  SAR 
theory. 

6.3.8  Comments  On  Bridge  Detection  and  Classification 

Reterence:  Figure  6.1.4 

The  San  Luis  Reservoir  SAR  photographs  are  very  good  tor  SAR  simulation 
evaluation.  The  potential  targets  in  these  photographs  include  two  dams,  a 
highway  bridge,  a  pumping  station,  two  transmission  lines,  and  a  switching 
yard.  The  elevation  in  the  top  photograph  goes  from  less  than  180  feet  to 
over  1100  feet.  In  the  top  photograph  the  O'Neill  Dam  is  in  the  top  ot  the 
photograph,  and  the  San  Luis  Dam  is  about  three-eighths  of  the  way  up  in  the 
photograph.  The  top  dark  area  is  the  O'Neill  Forebay.  The  bright  spots 
running  across  the  torebay  are  piers  that  contain  transmission  lines.  The 
white  dots  running  across  the  picture  are  transmission  line  towers.  The 
highway  bridge  is  just  a  dot  at  the  bottom  ot  the  forebay.  Note  the  water 
inlet  extending  out  from  the  San  Luis  Dam. 
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The  bottom  photograph  has  seven  times  the  resolution,  halt  the  range  and  twice 
the  depression  angle  ot  the  top  picture.  The  details  ot  the  highway  bridge 
are  apparent,  and  the  pumping  station  with  its  switching  yard  are  now  out  ot 
the  shadow  ot  the  San  Luis  Dam.  The  pumping  station  is  on  the  lower  lett  and 
the  switching  yards  are  to  the  right  ot  the  pumping  station.  The  transmission 
line  towers  are  not  as  apparent  as  in  the  top  photograph.  SAR  returns  are 
also  showing  up  trom  the  bottom  ot  the  reservoir.  The  two  photographs  seem  to 
be  processed  in  a  ditterent  manner  by  the  SAR. 

6.3.9  Comments  On  Haiwee  Reservoir 

Reterence:  Figure  6.1.5 

These  three  SAR  photographs  are  ot  the  northern  halt  ot  the  North  Haiwee 
Reservoir.  They  are  good  examples  ot  SAR  gain  setting  and  moving  targets. 
The  42- toot  resolution  picture  has  very  good  dynamic  range.  Interstate  395 
and  power  line  towers  run  north  and  south  at  the  left  ot  the  picture,  with  the 
highway  to  the  west  ot  the  power  line  towers.  Interstate  395  is  a  very  taint 
line,  but  the  motor  vehicles  and  buildings  are  its  signatures.  The  shadows 
and  water  are  quite  dark. 

The  17-toot  resolution  picture  has  a  high  gain  setting  and  the  shadows  have 
quite  a  bit  ot  noise.  The  bottom  ot  the  reservoir  shows  up  quite  well  and  the 
tace  ot  the  dam  gives  a  very  strong  radar  return.  The  power  line  towers  are 
almost  completely  eclipsed  by  the  high  gain  setting.  In  the  lower  lett  corner 
Interstate  395  can  be  seen  with  the  motor  vehicles  displaced  trom  the 
highway.  The  taces  ot  the  hills  give  very  strong  returns. 

The  8.5  toot  resolution  picture  has  a  lower  gain  setting,  but  the  bottom 
reservoir  is  still  visible.  The  tace  ot  the  dam  can  still  be  recognized  with 
the  reduced  return  and  shadow.  The  road  leading  to  the  base  ot  the  dam  is  now 
visible.  The  base  ot  the  dam  at  the  right  seems  to  have  some  man-made 
building.  There  are  other  bright  returns  that  can  only  be  explained  by  an 
on-site  inspection. 
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6.3.10  Comments  On  Airfield  Recognition 
Reference:  Figure  6.2.1 

Airports  are  ot  special  interest  to  the  military.  These  SAR  pictures  show  two 
airports.  The  China  Lake  Airport  is  modeled  in  the  Link- generated  data  base. 
The  other  airport  is  the  Inyokern  Airport.  Both  airports  are  in  an  isolated 
desert  area;  hence  they  are  very  prominent.  Both  pictures  show  saturation. 
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7.0  EMULATIONS  OF  SAR  IMAGES 

To  verity  that  today's  state  ot  the  art  in  electronics  is  adequate  to  build 
its  SAR  simulator.  Link  emulated  the  design  in  the  Link  Emulation  Laboratory. 
The  emulation  consisted  ot  a  FORTRAN  program  that  perlormed  the  same  function 
in  non- real  time  as  the  SAR  simulator  would  in  real  time.  The  generation  ot 
each  image  required  typically  20  to  40  minutes. 

Two  ditterent  types  ot  data  bases  were  used.  The  first  one  was  DMA  level  I, 
tor  the  area  ot  China  Lake  NAS,  California.  This  area  was  chosen  because  Link 
has  actual  SAR  inquiry  ot  the  area. 

In  addition  to  the  Level  I  data,  three  areas  covered  by  Level  X  were  also  used 
to  generate  SAR  images.  These  areas  were:  Scott  AFB ,  Lambert  Field,  St. 
Louis  Missouri,  and  Chester,  Illinois. 

As  anticipated,  DMA  Level  I  data  is  clearly  inadequate  tor  realistic  SAR 
simulation.  To  produce  realistic  SAR  images,  the  DMA  data  base  had  to  be 
enhanced  by  adding  buildings,  taxiways,  and  other  features  as  explained;  i.e., 
Section  6.2.1. 

In  this  section  the  emulated  SAR  images  will  be  shown  with  and  without  certain 
ettects.  They  clearly  demonstrate  that  synthetic  signature  generation,  the 
addition  ot  noise  and  various  other  ettects,  is  necessary  tor  realistic 
simulation. 

The  first  six  images  are  at  42  toot  resolution,  the  same  as  the  photograph 
made  by  the  F- 15  radar  tor  the  same  area,  as  shown  in  Figure  6.2-1.  The 
emulated  images  were  made  at  the  same  estimated  alrcratt  position,  depression 
angle  and  other  parameters  as  the  real  photograph.  A  short  explanation  ot 
each  image  follows. 

Figure  7.0-1  was  produced  trom  DMA  Level  1  data  by  using  present  radar 
landmass  simulation  techniques;  i.e.,  without  synthetic  signature  generation, 
noise,  overload,  etc.  Only  the  computational  accuracy  was  increased  to 
correspond  to  a  five  foot  least  significant  bit.  A  diagonal  line  in  the  upper 
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right  hand  corner  is  apparent  as  well  as  some  line  structure  in  the  lower 
right  hand  corner.  Both  are  due  to  problems  in  the  DMA  elevation  data  base. 
The  dark  line  is  due  to  an  "embankment"  being  created  by  the  elevation  data 
base.  Such  problems  are  generally  not  noticeable  with  real  beam  radars 
because  ot  the  lower  resolution.  The  single  large  teature  above  the  runway 
is,  in  tact,  a  collection  ot  buildings.  At  SAR  resolution,  the  individual 
building  would  be  distinguishable,  but  the  Level  I  data  base  does  not  provide 
such  detailed  inlormation. 

Figure  7.0-2  is  the  same  image  but  synthetic  signature  generation  (a  hardware 
subsystem  in  the  simulator)  has  been  added  to  break  up  the  ground  return  and 
the  return  trom  the  buildings. 

Figure  7.0-3  is  the  same  area  but  the  data  base  has  been  manually  enhanced  on 
the  basis  ot  other  data  as  explained  in  Section  6.2.1.  Figure  7.0-3  has  no 
synthetic  signature  generation,  while  Figure  7.0-4  uses  synthetic  signature 
generation. 

Figure  7.0-5  again  uses  the  unenhanced  data  base  but  with  synthetic  signature 
generation,  glitter  and  the  "ventetian  blind  ettect,"  caused  by  a  tall-ott  in 
the  antenna  gain  towards  the  edges  ot  the  swath.  1A  Venetian  blind  ettect  may 
also  be  caused  by  velocity  error  or  tootprint  error  as  explained  in  the 
report.)  Figure  7.0-5  is  an  enhanced  data  base  version  ot  Figure  7.0-6. 

The  next  tour  images,  Figures  7.0-7  through  7.0-10  are  at  a  resolution  ot  20 

teet  and  all  use  the  enhanced  data  base.  Figure  7.0-7  has  only  synthetic 

signature  generation.  Figure  7.0-8  has  receiver  noise  added.  Figure  7.0-9 
contains  the  ettect  ot  t liter  broadening  due  to  a  strong  return.  This  ettect 
occurs  in  both  the  azimuth  and  range  directions.  In  Figure  7.0-10 

noise/glitter  has  been  added,  while  Figure  7.0-11  has  t liter  broadening  and 

analog-to-digltal  overload. 

The  remaining  tour  images,  7.0-12  through  7.0-15  have  a  10  toot  resolution. 
The  enhanced  data  base  has  12  F/A-18  aircratt  lined  up  in  a  row,  parked 
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Figure  7.0-4  SAKE  AS  FIGURE  7.0-3  BUT  WITH  SYNTHETIC  SIGNATURE 

GENERATION  ADDED 
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Figure  7.0-5  SAME  AS  FIGURE  7.0-1  WITH  VARIOUS  SAR  EFFECTS 
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Figure  7.0-8  ENHANCED  DATA  IASI  VERSION  OP  FIGURE  7.0-5 
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Figur*  7.0-7  MMUALLY 
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Figure  7.0-11  SM  IMAGE  WITH  ADDITIONAL  SAR  EFFECTS 
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Figure  7.0-13  SM  IMAGE  WITH  NOISB 
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parallel  to  a  runway  (Figure  7.0-12).  All  images  have  synthetic  signature 
generation.  Figure  7.0-13  has  noise  added,  Figure  7.0-14  has  tiller 
broadening,  while  Figure  7.0-15  incorporates  the  ettects  ot  synthetic 
signature  generation,  filter  broadening,  and  noise/glitter. 

The  next  tour  images  (7.0-16  through  7.0-19),  were  generated  from  DMA  Level  x 
data.  The  tirst  ot  these.  Figure  7.0-16  is  a  false  color  map  ot  Scott  AFB. 
The  different  colors  correspond  to  different  surface  material  .  The  shades 
are  not  represent ive  ot  the  reflectance  values  magnitude.  Figure  7.0-17  shows 
a  radar  image  ot  a  portion  ot  Scott  AFB  at  a  resolution  ot  10  ft. 

Figure  7.0-18  is  the  SAR  image  of  Lambert  Field  at  a  resolution  ot  20  tt., 
while  Figure  7.0-19  shows  a  bridge  near  Chester,  Illinois. 
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8.0  CONCLUSIONS  AND  RECOMMENDATIONS 

On  the  basis  ot  this  SAR  Simulation  Study  and  other  extensive  research  on  the 

subject.  Link  has  reached  the  following  conclusions  and  recommendations. 

1)  Currently  existing  digital  landmass  simulators  are  inadequate  tor 
ettective  simulation  ot  synthetic  aperture  radar.  SAR  simulation 
requires  higher  accuracy  (at  least  a  live  toot  pixel  and  complex  SAR 
signature  generation,  neither  ot  which  are  available  in  present  systems. 
However,  the  present  state  ot  the  art  in  digital  technology  is  adequate 
tor  building  a  SAR  simulator. 

2)  Because  all  new  radar  systems  will  be  controlled  by  powertul  general 
purpose  computers  and  because  mode  selection  will  be  accomplished  through 
these  computers  via  programming,  it  is  very  likely  that  over  the  lifetime 
ot  the  radar  system  a  great  number  ot  significant  operational  changes 
will  take  place.  In  order  to  be  cost-ettective,  the  simulator  will  need 
to  have  a  similar  kind  ot  organization;  that  is,  to  primarily  be  a 
programmable  device  rather  than  a  hard-wired  one.  Only  system  operations 
that  are  part  ot  radar  physics  and  remain  invariant,  should  be  hard  wired 
in  the  simulator. 

3)  The  greatest  problem  in  synthetic  aperture  radar  simulation  is  due  to 
inadequacies  of  data  bases.  In  order  to  realistically  simulate  an  area 
to  ground  truth  accuracy,  a  Level  X  data  base  must  exist.  Because  ot  the 
very  high  cost  ot  mapping  an  area  at  Level  X  detail,  it  is  extremely 
unlikely  that  large  areas  will  be  available  at  this  level  ot  detail. 

In  order  to  achieve  high  fidelity  simulation,  Link  recommends  tour 
approaches  to  alleviate  the  data  base  problem: 

a.  When  Level  X  detail  is  available  in  the  vicinity  ot  navigation 
checkpoints  and  target  areas,  this  level  should  be  used. 
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b.  It  Level  X  detail  is  not  available  and  it  is  vitally  important  to 
have  actual  ground  truth,  manual  enhancement  ot  the  data  base  should 
be  used,  as  demonstrated  by  the  photographs  ot  the  China  Lake  area 
(reler  to  Figures  7.0-1  through  7.0-15).  It  ground  truth  is  not  ot 
great  importance.  Level  I,  II,  or  V  data  should  be  used,  combined 
with  synthetic  signature  generation  that  portrays  radar  signature  ot 
typical  teatures. 

c.  Even  with  SAR  signature  generation,  it  would  be  obvious  where  the 
boundary  ot  a  high  resolution  area  is  trom  the  roads,  streams,  and 
other  teatures  which  terminate  at  the  boundary.  To  make  this 
boundary  indiscernible,  Link  proposes  that  the  data  base  outside  the 
Level  X  area  be  augmented  by  artilicial  data.  This  artificial  data 
would  extend  roads,  streams,  and  selected  teatures  trom  the  Level  X 
area  by  using  a  random  data  base  generation  process  similar  to  the 
one  currently  employed  by  Link  in  visual  programs.  The  statistical 
characteristics  ot  the  artificial  data  would  correspond  to 
statistical  data  derived  trom  the  actual  area. 

d.  The  much  higher  resolution  ot  SAR  has  shown  anomalies  in  the  Detense 
Mapping  Agency  terrain  elevation  data  base.  These  anomalies  may  be 
due  to  problems  in  the  algorithms  used  in  the  conversion  ot  contour 
map  data  to  grid  elevation  posts,  or  merely  due  to  the  tact  that  the 
elevation  data  has  a  least  significant  bit  ot  one  meter.  This  causes 
quantization  ettects  which  result  in  "embankments"  or  the  appearance 
ot  terracing  ettects.  A  solution  must  be  found  either  by  the  Defense 
Mapping  Agency  or  by  the  simulation  manufacturer  to  eliminate  these 
spurious  ettects. 


On  the  basis  ot  studying  radar  photographs,  Link  concludes  that  a  number  ot 
special  ettects  and  malfunctions  should  be  simulated  tor  realistic  simulation 
ot  synthetic  aperture  radar.  The  following  items  should  be  included: 
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Azimuth  and  Range  Direction  Layover  -  Layover  in  both  azimuth 
and  range  direction  should  be  simulated.  (Range  direction 
layover  applies  equally  to  a  real  beam  radar  simulation.) 

Moving  Object  Displacement  -  Moving  objects  on  the  ground  are 
displaced  in  the  azimuth  direction  due  to  the  basic  assumption 
in  synthetic  aperture  radar  that  all  ground  objects  are 
stationary.  Simulation  should  include  displacement  ot  the 
object  and  the  correct  placement  ot  the  shadow. 

Velocity  Errors  -  Velocity  errors  must  be  simulated  because  the 
pilot  should  be  trained  in  recognizing  such  errors  and  taking 
corrective  action  by  recalibrating  the  inertial  navigation 
system.  A  small  velocity  error  would  result  in  feature 
displacement  in  the  azimuth  direction.  Large  velocity  errors 
can  completely  shitt  the  ground  map  outside  the  area  to  be 
viewed.  It  this  happens,  the  beam  width  ot  the  radar  may  be 
insufficient  to  illuminate  the  desired  area.  This  is  called 
footprint  error. 

Incorrect  Elevation  -  Simulation  ot  incorrect  assumed  elevation 
is  needed  because  the  radar  assumes  that  all  objects  on  the 
ground  are  in  a  plane  and  the  distance  between  the  plane  and 
the  aircraft  is  known.  Should  these  assumptions  be  untrue,  the 
map  will  be  distorted  by  effects  similar  to  velocity  errors. 

Motion  Compensation  or  Acceleration  -  Error  due  the  flexing  ot 
the  aircraft  or  the  aircraft  not  flying  in  a  straight  line 
implies  that  the  velocity  at  the  inertial  navigation  system  may 
not  be  identical  to  the  velocity  at  the  radar  antenna. 
Depending  on  the  particular  radar,  various  methods  may  be  used 
to  ~orrect  this  disparity.  The  failure  ot  this  correction 
should  be  simulated. 
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(6)  Antenna  Side  Lobe  Etlect  -  The  antenna  side  lobe  ettect  may 
produce  an  unwanted  input  with  synthetic  aperture  radar. 
Features  withih  this  side  lobe  will  be  superimposed  on  the 
regular  image.  Also,  in  the  range  direction,  the  Barker  code 
will  produce  an  ettect  similar  to  that  of  the  antenna  side 
lobes.  Typically,  a  strong  signal  will  recur  in  the  range 
direction  at  fixed  intervals  trom  the  actual  location  ot  the 
return.  The  pilot  should  be  trained  to  recognize  the 
appearance  ot  these  ettects. 

(7)  Filter  Broadening  -  Strong  returns  have  a  tendency  to  saturate 
the  digital  t liter  and  certain  teatures  may  appear  to  be  larger 
than  their  actual  size.  Similarly,  filter  saturation  may 
result  in  the  repetition  ot  a  strong  signal  in  the  range 
direction. 

(8)  Analog- to-Digital  Overload  -  An  overload  in  the 
analog- to-digital  converter  will  produce  harmonics.  These 
harmonics  manitest  themselves  as  spurious  returns  on  an  image. 
The  simulation  ot  analog- to-digital  overload  should  be 
implemented  whenever  such  an  overload  can  occur. 

(9)  Glitter  -  The  well  known  glitter  ot  synthetic  aperture  radar 
can  be  caused  by  the  coherent  properties  ot  the  transmission 
and  by  objects  on  the  ground  (such  as  buildings)  that  act  as 
corner  retlectors.  Glitter  should  be  simulated  to  maintain 
high  fidelity  simulation  so  that  pilots  would  not  be  trained  on 
an  idealized  image.  An  idealized  image  (i.e.,  no  glitter) 
could  give  the  pilot  a  false  impression  of  the  correct 
appearance  ot  the  image  and  effect  transfer  of  training. 

(10)  Multilook  -  With  most  synthetic  aperture  radar  systems,  several 
ground  map  returns  are  averaged  under  the  control  ot  the  pilot 
to  reduce  glitter  and  noise.  Unless  the  display  equipment  is 
actual  aircraft  hardware,  the  simulator  must  perform  this 
multilook  averaging  function. 
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(11)  Weather  Simulation  -  Depending  on  the  wavelength  ot  the  radar, 
clouds  or  moisture  will  generate  strong  returns  and  also  cause 
absorption.  Absorption  will  cause  a  loss  ot  signal  strength 
and  a  shadow  in  the  area  behind  the  cloud.  Weather  simulation 
may  be  necessary  depending  on  the  wavelength  ot  the  radar. 
Moving  clouds  should  be  displaced  in  a  manner  similar  to  moving 
objects  on  the  ground. 

The  above  list  ot  ettects  are  those  that  should  be  simulated  tor  most 
synthetic  aperture  radars.  Depending  on  the  particular  radar  or  the  mission 
ot  the  aircraft,  all  ot  these  ettects  may  not  be  required  tor  ettective 
training.  Since  SAR  simulation  is  a  new  technology  and  consequently  the 
requirements  tor  a  SAR  simulator  are  full  ot  uncertainties,  Link  recommends 
that  the  government  procure  a  research  SAR  simulator. 

The  purpose  ot  this  simulator  would  be  to  develop  methods  tor  advanced  senior 
simulation  and  to  tully  explore  the  specifications  ot  new  simulators.  A 
secondary  purpose  ot  this  research  simulator  would  be  to  aid  in  the 
development  ot  air-to-ground  tactics. 


I 
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DEFINITION  AND  SYMBOLS 

Ff  Transmitted  frequency 

Fg  Aperture  function 

G'c  Gain  function  (due  to  gain  control,  weather,  and 
atmospheric  attenuation) 

Receiver  antenna  gain 
Transmitter  antenna  gain 

//  Height 

i  Imaginary  constant  \T  1  "  (physics  usage) 

j  Imaginary  constant  \J~  T  “  (electronic  usage) 

K  Selectable  constant 

k  Wave  number 

Kf  Aperture  modification  factor 

hf  Earth  radius  modification  factor 

Aa  Selectable  constant  for  SAR  resolution 

KF  Phase  shift  limit 

hH  cos(6w)  cos(0w) 

L/p  Isodoppler  line 

O  Orientation 

P  Phase  shift 

PfK  Pattern  factor  along  the  receiving  path 
PfT  Pattern  factor  along  the  transmitting  path 
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DEFINITION  AND  SYMBOLS 

Pff  Received  power 

Pf  Transmitted  power 

R  Range  for  radar 

R  Distance  from  aperture  to  image  plane  (physic 

usage) 

Rdb  Data  base  resolution 

Rds  Aircraft  display  resolution 

R e  Earth  radius 

Rec  Modified  earth  radius  for  atmospheric  refraction 

Rfu  Range  forward  mapping  due  to  elevation 

RB  Range  rate 

Rsaf  Synthetic  aperture  radar  resolution 

Rsl  Slant  range 

Rsum)  Range  coordinate 

Rtl  Range  at  time  t, 

Rt 2  Range  at  time  <2 

Rw  Range  walk 

/?,  Range  at  position  one 

R2  Range  at  position  two 

Rues  Resolution  in  length 

Rdg  Delta  ground  in  length 

Rslt  Slant  range  top  of  target 

sin  Sine 


Samt  Received  signal  amplitude 
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DEFINITION  AND  SYMBOLS 


/shadow  effect 

S/p  In  phase  component  of  received  signal 
Sq  Quadrature  component  of  received  signal 

Sr  Received  signal 

Spp  Received  signal  in  frequency 

ST  Transmitted  signal 

Intensity  spread  function 

t  Time 

tan  Tangent 

Tcs  Target  cross  section 

(d  Delta  time 

Tret  Target  reflectance  coefficient 

u  Dimension  in  the  aperture 

«(,)  Complex  amplitude  distribution  of  a  point 

Fyic  Aircraft  velocity 

Vc  Closing  velocity 

Velocity  error 

Vgac  Velocity  from  aircraft  to  ground 

Vgt  Velocity  of  ground  target 

Vp  Velocity  between  the  aircraft  nadir  point  and  a 
point  on  the  ground 
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